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INTRODUCTION 

Although  the  Nyquist  criterion  [1]  has  been  known  for  over  half  a 
century,  it  has  resisted  generalization  until  recently.  Interestingly, 
those  generalizations  which  have  been  formulated  retain  the  simple 
graphical  character  of  the  classical  test,  even  when  one  is  studying 
systems  defined  on  abstract  spaces.  The  earliest  generalizations  of 
the  Nyquist  criterion  were  the  Circle  and  Popov  criteria  formulated 
in  the  early  sixties  as  nonlinear  and/or  time-variable  perturbations 
of  the  classical  test  [2],  [3].  More  recently  MacFarlane  [4],  and' 

Barman  and  Katznelson  [5]  have  extended  the  test  to  the  case  of  fre- 
quency response  matrices  while  one  of  the  authors  has  formulated  a 

Nyqui st-like  sufficiency  condition  for  Lipschitz  continuous  operators 
on  abstract  spaces  [6].  Finally,  in  a recent  paper  the  authors  gave  a 
stability  test  for  multivariable  digital  filters  which  was  formulated 
in  terms  of  a continuum  of  Nyquist  plots  [7],  [8].  In  all  cases  the 
tests  remain  simple  graphical  conditions  on  the  complex  plane.  The 
resultant  criteria  yield  necessary  and  sufficient  conditions  in  the  case 
of  linear  time-invariant  systems  (including  the  multivariable  and  matrix 
generalizations)  and  sufficient  conditions  in  the  cases  of  nonlinear 
and  time-variable  systems. 

The  thrust  of  this  paper  is  to  show  that  with  a slight  modification, 
the  continuum  of  Nyquist  plots  used  in  the  multivariable  Nyquist  test 
of  [7]  can  be  reduced  to  a classical  single  variable  Nyquist  plot  plus 
a test  to  verify  that  the  filter  has  no  poles  in  the  region  |zj|  3 |z2| 
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= ...  3 I zn | = 1,  the  multivariable  analog  of  the  iu-  axis. 

Although  every  attempt  is  made  to  minimize  details  the  theorem 
illustrates  the  essential  algebraic  topological  nature  of  the  Nyquist 
criterion,  with  homotopic  arguments  playing  a significant  role  in 
derivation.  For  a more  detailed  discussion  of  the  algebraic  topological 
nature  of  Nyquist  theory  the  reader  is  refered  to  [7],  [8]  and  [9]. 

In  particular,  [9]  gives  an  algebraic  topological  derivation  of  the 
classical  Nyquist  criterion. 

In  the  following  section  several  Hurwitz-like  stability  tests 
are  reviewed  and  a new  test  of  the  Hurwitz-type  test  is  formulated. 

In  the  third  section  a homotopic  interpretation  of  the  classical  Nyquist 
criterion  is  formulated,  this  being  used  to  derive  the  desired  multi- 
variable  Nyquist  criterion  from  the  Hurwitz  conditions  in  the  fourth 
and  fifth  sections.  Finally  some  examples  of  the  theory  are  given  in 
section  6. 


Hurwitz-like  Tests 

Denote  the  vector  space  of  complex  n-tuples  by  ?n.  For  the 
purposes  of  our  multivariable  stability  theory,  there  are  five  inter- 
esting subsets  of  0n.  First  there  is  the  polydisk,  defined  as 

Pn  = {(z1,...,zn)  in  ln  1 z - ( <1,  i = 1 , . . . ,n}  (1 ) 

It  plays  the  same  role  in  the  multivariable  theory  as  the  unit  disk 
(or  right  half  plane)  in  single  variable  theory. 

Next,  there  are  three  separate  notions  of  the  boundary  of  Pn. 

All  are  necessary  for  the  theory  of  this  paper.  First  is  the  dis- 
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tinguished  boundary  denoted  by  Tn  where 

Tn  - {(zt,...,zn)  in  ?n  [ |z1 | =1,1=  l,...,n}  (2) 

Tn  serves  as  the  multidimensional  analog  of  the  iu-axis.  In  particu- 
lar the  frequency  response  [10]  of  a digital  filter  is  the  evaluation 
of  its  transfer  function  over  Tn.  Geometrically  Tn  is  an  n-dimensional 
torus  which  reduces  to  the  unit  circle  of  the  complex  plane  in  the 
single  variable  case.  Tn  is  a "boundary"  for  Pn  in  the  sense  that  it 
is  a subset  of  Pn  for  which  all  coordinants  of  Pn  simultaneously  take 


on  extremal  values. 


A second  notion  of  boundary  for  Pn  is  defined  by  requiring  only 
that  n-1  coordinats  take  on  extremal  values.  This  boundary  set  is 
denoted  by  Mn  where 

Mn  = {(z1 zn)  in  f |}Zi|  = 1,  i = l,...,k-l,k+l, 

• • • » n » | zk I <_  1 } (3) 

and  where  k ranges  from  1 through  n. 

The  final  notion  of  boundary  requires  that  at  least  one  of  the 
coordinates  take  on  extremal  values.  This  notion  is  the  usual  topo- 
logical boundary  since  it  coincides  with  the  usual  concept  of  boundary 
of  the  set  Pn  in  the  sense  of  point  set  topology  [11].  The  topological 
boundary  is  denoted  by  Bn  where 


B = {(zlf...,zn)  in  l | | zi  [ il,  i = 1 ,...,n  and  |zk|  = 1 for 

some  k)  (4) 

Finally  we  define  a subset,  Hn,  of  Pn  whose  relevence  to  the 
stability  problem  was  originally  indicated  by  Huang  [12]. 
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Hn  = {(z1>...,zn)  in  £n  |Zi|  = 1,  1-1 k-1 ; [zk|  <1;  z t = 


0,i=k+l ,. ..  ,n)  f (5) 

Here  k varies  from  1 through  n.  Note,  that  M1  = H1  = P1  and  T1  = B1 

hence  these  sets  become  redundant  and  all  reduce  to  either  the  unit 

disk  or  unit  circle  in  the  single  variable  case. 

Shanks  [13]  was  the  first  to  give  a Hurwitz-like  test  for  the 

stability  of  multidimensional  digital  filters.  His  condition  states 

essentially  that  the  filter  transfer  function  must  have  no  poles  in  Pn. 

In  the  single  variable  case,  the  pole  set  of  a transfer  function  \s 

discrete.  However,  in  the  case  of  higher  dimensional  filters,  the 

pole  set  is  an  infinite  continuum.  Using  this  fact,  Huang  [12]  showed 

2 

that  a transfer  function  has  a pole  in  P if  and  only  if  it  also  has  a 
2 

pole  in  H . This  is  not  toimply  that  the  only  poles  of  the  transfer 
2 

function  lie  in  H but  rather  that  the  pole  set  is  so  large  that  it  cannot 

2 2 
pass  through  P without  intersecting  the  subset  H . Anderson  and  Jury 

[14]  extended  Haung’s  theorem  to  the  n-dimensional  case  by  showing 

that  a transfer  function  has  a pole  in  Pn  if  and  only  if  it  has  a pole 

in  Hn.  The  proof  of  Haung's  theorem  and  its  generalization  is  tediously 

straightforward  but  requires  a clever  application  of  the  maximum  modulus 

theorem. 

A result  somewhat  similar  to  Huang's  can  be  formulated  in  terms 
of  the  topological  boundary.  To  derive  such  a condition,  one  exploits 
the  fact  that  the  pole  set  of  a multivariable  (n  ^2)  rational  function 
is  an  infinite  continuum  (more  precisely  no  connected  component  of  the 
pole  set  is  compact  [15]).  As  such,  the  only  way  the  pole  set  can 
intersect  Pn  is  if  it  crosses  the  topological  boundary,  Bn.  This 


O 


Implies  that  a transfer  function  has  a pole  set  intersecting  Pn  if 
and  only  if  the  pole  set  has  a non-void  intersection  with  Bn. 

Now  observe  that  Bn  can  be  viewed  as  the  union  of  a family  of  (n-1) 
variable  polydisks  (parameterized  by  k and  the  value  of  zk,  |zk|  =1). 
Hence  the  above  argument  can  be  repeted  to  show  that  the  transfer 
function  has  a pole  set  intersecting  the  topological  boundary  of  such 
an  (n-l)-variable  polydisk  if  it  has  a pole  set  intersection  Pn. 

Upon  iterating  the  argument  (n-1)  times  and  eliminating  redundant  sets, 
one  eventually  arrives  at  the  following  condition:  the  transfer 
function  has  a pole  in  Pn  if  and  only  if  it  has  a pole  in  Mn. 

The  above  various  Hurwitz-like  stability  tests  for  multivariable 
digital  filters  are  summarized  as  follows: 

Therorem  1:  Let  a causal  multidimensional  digital  filter  be  character- 
ized by  a rational  transfer  function  in  several  complex  variables. 

Assume  the  numerator  and  denominator  polynomials  are  relatively  prime. 
Then  the  following  are  equivalent  stability  conditions: 

i)  the  pole  set  of  the  transfer  function  has  a null 
intersection  with  Pn. 

ii)  the  pole  set  of  the  transfer  function  has  a null 
intersection  with  Bn. 

iii)  the  pole  set  of  the  transfer  function  has  a null 
intersection  with  Hn. 

iv)  the  pole  set  of  the  transfer  function  has  a null 
intersection  with  Mn. 

The  easiest  way  to  evaluate  the  stability  tests  based  on  the 
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above  conditions  is  by  a comparison  of  the  (topological ) dimension 
of  the  sets  where  one  checks  for  the  existence  of  poles.  In  particular, 

Pn  is  2n-dimensional , Bn  is  (2n-l )-dimensional , while  Hn  and  Mn  are 
both  (n+1 )-dimensional . Again  realize  that  the  equivalences  of  Theorem 
1 follow  from  the  fact  that  the  pole  set  is  an  infinite  continuum 
whose  complex  dimension  is  (m-1)  where  m is  the  number  of  complex 
variables  of  the  specific  function.  Finally  observe  that  a pole  is 
implicitly  used  to  mean  a specific  point  in  the  "pole  set." 

Nyquist  Theory 

The  task  of  this  section  is  to  construct  the  concepts  of  a 
Nyquist  contour,  a Nyquist  plot,  encirclement,  and  degree,  all  in 
topological  terms  which  are  thus  extendable  to  the  multivariable  case. 

All  of  our  conditions  will  be  stated  in  terms  of  the  zero  set  of  a 
relatively  prime  denominator  polynomial  of  a transfer  function.  Hence 
we  will  deal  exclusively  with  polynomials  in  several  complex  variables 
rather  than  rational  transfer  functions. 

Traditionally  engineers  view  the  Nyquist  contour  as  a subset  of 
the  complex  plane.  This  point  of  view  is  somewhat  erroneous.  Mathe- 

I; 

matically  speaking  the  Nyquist  contour  (the  usual  closed  semicircle, 
the  imaginary  axis,  or  the  unit  circle)  is  a continuous  map  (of  bounded 
variation)  from  T1  to  £n.  The  image  of  this  map,  called  the  trace  of 
the  map,  is  the  traditional  engineering  notion  of  the  Nyquist  contour. 

In  this  paper  T1  is  the  unit  circle  of  the  complex  plane.  In  the  single 
variable  case  (classical  digital  filter  stability)  one  works  with  a 
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"Nyquist  contour"  defined  hy  r(a)  = a = exp(ie)  for  0 <_  0 < 2*,  0 =■  arg(a). 
In  the  multivariable  case,  the  map  will  take  on  values  in  £n  forcing  the 
"Nyquist  contour"  to  be  a more  involved  entity.  Observe  that  we  are  tak- 
ing liberties  with  the  calssical  definition  of  the  Nyquist  contour  and 
plot.  In  an  abstract  sense,  there  is  no  essential  difference  although 
the  specific  applications  (classical  feedback  stability  or  presently 
digital  filter  stability)  are  somewhat  alien. 

In  this  paper  a Nyquist  plot  is  defined  as  the  composition  of  the 
Nyquist  contour,  r,  with  a polynomial  in  several  complex  variables, 
r,  as  per  Figure  1.  Note  that 


Figure  1:  The  Nyquist  plot  as  a composition  of  maps. 

the  polynomial  map  from  £n  to  l is  analytic  and  it  is  this  property  which 
makes  the  theory  go.  Thus  the  Nyquist  plot  is  a continuous  map  of  bounded 
variation  from  T^  to 

The  concept  of  encirclement  is  intimately  related  to  the  topological 
concept  of  homotopy.  Keeping  this  association  in  mind,  let  X be  an 
arbitrary  topological  space  and  let  u and  y be  continuous  ^unctions 


pm 


of  bounded  variation  defined  on  T1  with,  values  in  X:  y.T^X  and  xtT^X. 
The  maps  y and  x are  said  to  be  homo topic  if  there  exists  a continuous 
map,  0,  defined  on  the  product  space  t’xI,  I * [0,7],  with  values  in  X 
such  that  0(a,O)  = y(a)  and  0(a,l)  = x(a).  In  essence,  a homotopy  is 
a continuous  deformation  of  the  curve  y into  the  curve  X.  This  concept 
defines  an  equivalence  relation  on  the  set  of  continuous  maps  from  T1 
to  X— i.e.  two  curves  are  equivalent  if  one  can  be  continuously  deformed 
into  the  other.  A curve  is  said  to  be  homotopically  trivial  if  it 
is  homotopic  to  a constant  map.  Note  that  the  use  of  an  abstract 
topological  space,  X,  in  the  definition  of  homotopy,  is  fundamental  'to 
the  concept,  since  all  curves  with  values  in  0n  or  Rn  are  homotopically 
trivial.  Although  we  are  interested  in  the  properties  of  functions 
defined  on  0n,  a number  of  non-trivial  topological  spaces  arise  in  our 
analysis.  In  particular  the  torus,  Tn,  and  the  punctured  plane.  Here  th 
concept  of  encirclement  may  be  defined  for  maps  taking  their  values  in 
0 - (0).  This  is  a highly  nontrivial  space  in  which  the  distinct 
equivalence  classes  of  homotopic  maps  can  be  indexed  by  the  integers 
corresponding  to  the  number  of  times  a curve  encircles  the  point  zero. 
This  number  is  termed  the  degree,  n(y,o),  of  the  map,  y,  and  can  be 
computed  by  the  formula  [11] 

n(y,0)  = 2^  / (z-0)-1dz  (6) 

Y 

where  y is  the  curve  in  question  and  y does  not  take  on  the  value  zero. 
The  desire  is  to  make  a binary  decision  on  whether  or  not  the  map  en- 
circles zero.  This  may  be  defined  in  purely  homotopic  terms  by  saying 
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that  a map  u:T^-*£  - {0}  does  not  encircle  zero  if  it  is  homotopically 
trivial  (homotopic  to  a constant  map).  The  concept  may  then  be  extend- 
ed to  maps  taking  values  in  g which  do  not  pass  through  the  point  zero 
by  viewing  them  as  maps  with  values  in  £ - {0}. 

The  significance  of  these  homotopic  concepts  in  algebraic  topology 

[11]  is  due  to  the  fact  that  the  equivalence  classes  of  homotopic  maps 

form  a group,  *(X),  where  X is  the  space  in  question.  ir(X)  is  called  the 

fundamental  group  of  X.  More  precisely,  if  we  have  two  maps  piT1  -*■  X 

and  x-.T1  X,  such  thaty(l)  = x(l),  their  concatenation  y * x:f  -»•  X 

is  defined  by  • 

r 2 

V X(a  ) 0<  arg(a)  < tt 

[y*x](a)  = -j 

[_u(a  ) ir  <arg(a)  < 2u 

Intuitively,  u*x,  is  a curve  which  first  follows  x as  arg(a)  goes  from 
0 to  tt  and  then  follows  y as  arg(a)  goes  from  tt  to  2tt.  Moreover,  since 
u(l)  = x(l),  then  y*x  is  continuous  if  u and  X are  continuous.  Clearly 
concatenation  is  invariant  under  homotopic  equivalence  [11].  (i.e.  if 

u-j  is  homotopic  to  U2  an<^  1S  homotopic  to  t^ien  1S  homo- 
topic to  V2  * As  suc^  ^ defines  a binary  operation  on  the  equivalence 

classes  of  homotopic  curves  jkT1  -»■  X with  a fixed  value  for  y(l).  Thus  the 
group  operation  of  tt(X)  is  concatenation  [11].  In  classical  algebraic  topo- 
logy, the  properties  of  ir(£n)  = 0 since  all  maps  taking  values  in 
are  homotopically  trivial.  On  the  other  hand  n(£  - {0})  is  isomorphic 
to  the  additive  group  of  integers  where  the  "degree  function"  is  the 
isomorphism. 


t 
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Using  this  machinery  we  now  formalize  a statement  of  the  "classical " 
Nyquist  criterion.  Consider  the  "obvious"  Nyquist  contour,  r:T*  -*■  l 
defined  as  r(a)  = a = exp[i  arg(a)]  with  the  corresponding  Nyquist  plot 
A°r:T1  -»■  J5  where  A is  a polynomial  in  one  variable. 

Theorem  2:  (Nyquist)  Let  A be  a polynomial  on  0.  Then  A has  no  zeros 
in  P1  (the  unit  closed  disk)  if  and  only  if  A°r  does  not  pass  through 


nor  encircle  zero. 


Main  Theorems 


Here  the  multivariable  Nyquist  theory  is  derived  from  condition 
(iv)  of  Theorem  1.  Let  a causal  digital  filter  transfer  function  be 
H(z-j»...  ,zn)  = B(z-j ,. ..  ,zn)/A(z-j ». . . ,zn)  where  A and  B are  relatively 
prime.  The  system  so  characterized  is  stable  if  and  only  if  A has  no 
zeros  in  Mn— i.e.  the  zero  set  of  A does  not  intersect  Mn.  Now  Mn 
can  be  expressed  as  a union  of  single  variable  polydisks  as  follows. 
First,  for  any  given  set  of  (n-1 ) elements,  c^,  of  T1 , indexed  by  the 
integers  l,2,...,k-l,  k+l,...,n,  embed  a single  variable  polydisk  into 
£n  as 

P1(a1,...,ak_1,*,ak+1,...,an) 

{(Zi ....  ,zn)  in  £ |zi  - a.j , i = l,...,k-l, 

k+1 n;  | zk|  <_  1)  ' (8) 

By  comparison  with  equation  3,  one  can  verify  that 
Mn  = U ot,  U pVi 0k-V‘»ak+r*“*an) 

t i l » 


Hence  the  digital  filter  is  stable  if  and  only  if  A has  no  zeros  in 
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each  of  the  polydisks  Moreover,  since 

these  are  dependent  on  only  one  coordinant,  one  may  test  for  zeros  of 
A in  (a-j ,•  *0^1 . - ,an)  by  sequentially  testing  for  zeros 
of  the  single  variable  polynomial  A(a-j .... » . . . >an)  in  the 
•single  variable  polydisk  P1  as  defined  in  the  introduction.  Each  such 
test  can  be  executed  using  the  "Nyquist  theorem."  This  yields  a stability- 
test  for  a function  of  several  variables  which  takes  the  form  of  a con- 
tinuum of  classical  Nyquist  plots. 

Lemma  3:  Let  A be  a polynomial  mapping  £n  to  £n.  Then  A has  no  zeros 
in  Pn  if  and  only  if  each  of  the  Nyquist  plots  for  the  family  of  single 
variable  polynomials. 


A(a-, , . . . jctk_-j  ,zk>  ak+l  **  * • ,cxn^  ^or  and  c-j  in  T^ , do  not 

pass  through  nor  encircle  zero. 

Lemma  3 is  essentially  equivalent  to  the  condition  formulated  in 
references  [7]  and  [8].  However,  the  present  condition  arose  from  the 
fact  that  A has  no  zeros  in  Mn  whereas  the  previous  test  grew  from  the 
fact  that  A has  no  zeros  in  Hn.  As  was  shown  in  references  [7]  and  [8] 
the  lemma  can  be  implemented  as  a practical  test  in  the  two  dimensional 
case.  Here  one  simply  chooses  a finite  set  of  a's  in  T1  and  plots  the 
corresponding  Nyquist  loci.  Since  T1  is  compact,  this  discretization 
can  be  made  to  yield  as  much  accuracy  as  desired.  Unfortunately,  in 
the  multivariable  case  our  family  of  Nyquist  plots  is  parameterized 
by  "n"  (n-1 )-dimensional  tori. 

The  purpose  of  the  first  main  theorem  is  to  show  that  the  family 
of  Nyquist  plots  of  Lemma  3 is  reducible  to  "n"  classical  single  variable 
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Nyquist  plots  and  one  further  test.  To  this  end  we  reformulate  the 
classical  Nyquist  contour  used  to  test  the  single  variable  function 
A(<^  , . . . »ctjc_|  »2)c>a|c+i  »•  • • *0^)  as  a Nyquist  contour  taking  its  values 
in  0n;  define  the  Nyquist  contour,  rtc^j ..  ,ak _1  ,•  ,ak+1 . ,a 
by  the  equality 

r(c,l’-*-»ak-i»ot»ak+l an^ 

= (a1,...,ak_1,ct,ak+1 aj  (9) 

Clearly,  the  Nyquist  plot  Aor(ai , . . . ,ak_.j  ,a,ak+1 ,. . . ,an)  coincides 
with  the  Nyquists  plots  of  Lemma  3.  Thus  Lemma  3 can  be  reformulated 
in  terms  of  these  plots.  The  key  attribute  of  these  multivariable 
Nyquist  contours  is  that  the  entire  family  of  Nyquist  contours  for  a 
fixed  k are  homotopical ly  equivalent. 

Lemma  4:  For  any  given  set  of  a^s  in  T1 , i = 1 , . . . ,k-l  ,k+l . . ,n, 
the  Nyquist  contour  rC«i » . • • »“k_i »* »ak+] * • • • >“n)  1S  homotopic  in  Tn 
to  the  Nyquist  contour  r(l ,. . . ,1 ,• ,1 ,. . . ,1 ) . 

Proof:  Consider  the  homotopy  0:T^xI->Tn  defined  by 

0(a,t)  = (exp[i01(l-t)],...,exp[iek-1(l-t)],a,exp[i9k+1(l-t)], 

...,exp[i9n(l-t)])  (10) 

where  ^ - argCa^  for  any  set  xi  eT1 , i = 1 ,. . . ,k-l  ,k-M n.  Here 

0(a,t)  is  in  Tn  for  all  a and  t.  Moreover 

0(a,O)  = [exp(i  e-| ),..., exp(iek_-,), a, 

exp(i0k+1),...,expCi9n)] 

s r(ai  ♦•  ••  >ak-l  ,a»ak+r' **»an^ 
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and 

0(a,l)  = 0 l»a,l,...,l)  = rO l,a,l,...,l)  02) 

Hence  0 is  the  desired  homotopy. 

Theorem  3:  Let  A be  a polynomial  on  0n.  Then  A has' no  zeros  in  Pn 
if  and  only  if: 

(i)  A has  no  zeros  on  T*1,  and 

(ii)  The  Nyquist  plots  for  th'e  single  variable  functions 

A(1 ,. . . ,1 ,zk,l , . . . ,1 ) k = 1 n do  not  encircle  zero. 

. 

Proof:  Since  the  images  of  the  family  of  Nyquist  contours  defined 
by  equation  9 cover  Tn,  the  fact  that  none  of  the  Nyquist  plots  of 
Lemma  3 go  through  zero  implies  that  A has  no  zeros  on  Tn.  This 
verifies  the  necessity  of  condition  (i).  To  verify  the  necessity  of 
condition  (ii)  observe  that  the  "n"  Nyquist  plots  of  condition  (ii) 
are  a subset  of  the  family  of  Nyquist  plots  of  Lemma  3. 

To  verify  the  sufficiency  of  the  theorem,  observe  that,  if  A 
has  no  zeros  on  T , then  A restricted  to  Tn  is  a continuous  map  from 

' 

T to  0 - {0}.  Now,  since  the  continuous  images  of  homotopic  maps 
are  homotopic,  the  fact  that  the  Nyquist  contours  r (ai  * . - . »ctk 
ak+l  * " ’ *an ) anc*  >--»»l)  are  homotopic  for  any  fixed  k 

with  a.j  in  T » i = 1 . . ,k-l ,k+l implies  that  the  Nyquist 

plots  Aor(ajc_i>*,ak+i,...,an)  and  Aor(l ,. . .',1 , • ,1 , . . . ,1 ) are  homo- 
topic in  0 - {0};  hence  all  such  Nyquist  plots  encircle  zero  if  and 
only  if  the  Nyquist  plot  for  Aor(l,...,l,*,l,...,l)  encircles.  As 
such,  if  A has  no  zeros  on  Tn  we  are  assured  that  none  of  the  Nyquist 
plots  of  Lemma  3 go  through  zero,  whereas  if  the  ”n"  Nyquist  plots 
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Aor(1,...,l,»,l,...,l)  do  not  encircle  zero,  then  none  of  the  Nyquist 
plots  of  Lemma  3 encircle  zero.  With  the  final  observation  that  the 
Nyquist  plot  for  the  single  variable  function  A(l,...,l,z  ,1,...,1) 

coincides  with  Aor(l 1, .,1)  this  verifies  the  sufficiency 

of  the  theorem. 

Surely  this  theorem  is  a true  generalization  of  the  classical 
Nyquist  theorem  in  that  it  tests  for  stability,  using  only  distinguished 
boundary  (iu-axis)  information.  Moreover,  the  test  is  n-dimensional , 
hence  superior  to  any  of  the  Hurwitz-type  tests. 

» 

Intuitively  speaking  the  result  is  both  surprising  and  expected. 

It  is  surprising  because  one  tests  for  zeros  of  an  n-variable  function 
using  single  variable  Nyquist  plots  as  opposed  to  some  type  of  n-dimen- 
sional  encirclement.  It  was  expected,  however,  since  a polynomial  con- 
tains a finite  amount  of  information  (a  finite  number  of  coefficients), 
so  that  only  a finite  number  of  tests  need  be  executed.  In  this 
light  the  condition  of  Lemma  3 seemed  superfluous. 

Again,  Theorem  3 is  aesthetically  pleasing  since  it  uses  only 
frequency  response  information.  However,  by  cleverly  considering 
the  implications  of  this  information  as  per  [12],  [14],  one  may  concoct 
an  equivalent  test.  Essentially  the  test  will  be  a consequence  of 
condition  (iii)  of  Theorem  1 and  hopefully  will  be  easier  to  implement. 
Theorem  4:  Let  A be  a polynomial  mapping  £n  to  ?n.  Then  A has  no 
zeros  in  Pn  if  and  only  if 

(i)  A has  no  zeros  on  Tn,  and 

Cii)  The  Nyquist  plots  for  the  single  variable  functions 

A(1 ,...,1 jZ^, &,..., 0)  k = l,...,n,  do  not  encircle  zero. 
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Proof;  Note  first,  that  from  [7],  and  [8],  it  is  known  that  A has  no 
zeros  in  Pn  if  and  only  if  the  Nyquist  plots  of  A(a-|  ,...,0^  1 .z^.Q,.. . ,0) 
for  k = l,...,n  and  for  each  and  every  a.,-  in  do  not  encircle  nor  pass 
through. zero. 

Since  A has  no  zeros  in  Pn,  it  immediately  follows  that  A has 
no  zeros  on  Tn.  Moreover,  the  Nyquist  plots  of  condition  (ii)  are  a 
subset  of  the  Nyquist  plots  of  the  stability  result  noted  above.  Thus 
the  forward  direction  is  shown. 

Conversely  suppose  condition  (.i)  and  condition  (.ii)  hold.  Since 
A has  no  zeros  in  Tn,  Lemma  4 guarantees  that  the  Nyquist  plots  of* 
A(ap.. . ,an_i  ,zn)  for  each  and  every  in  are  homotopic  to  one 

another.  In  particular  they  are  homotopic  to  A(1 l,zn).  Thus 

any  member  of  this  family  of  Nyquist  plots  encircles  zero  if  and  only 
if  A(1 ,. . . ,1 ,zn)  encircles  zero. 

Now  if  the  Nyquist  plot  of  A(1 ,. . . ,1 ,zn)  does  not.  encircle  zero, 
then  A(o-j,...,an_-|,0)  f 0 for  each  and  every  in  . In  other  words 
A(z-j ,. . . ,zp_i  ,0)  has  no  zeros  on  the  n-1  dimensional  torus;  {(z^,...,z  ) 
eCn  1 1 | = 1 ‘»i  = l,...,n-l;zn  = 0}-  Repeating  the  above  arguments  we 

conclude  that  the  Nyquist  plot  of  A(1 ,. ..  ,1  ,zn_-j  ,0)  encircles  zero  if 
and  only  if  each  of  the  Nyquist  plots  of  A(a^ ,. . . ,an_2,an  ^,0)  encircles 
zero.  Continuing  in  this  fashion,  one  verifies  that  the  conditions  of 
the  theorem  are  equivalent  to  the  condition  of  references  [7],  [8]  as 
stated  at  the  beginning  of  the  proof.  The  statement  of  the  theorem 
now  fol 1 ows . 


Observe  that  this  theorem  does  not  explicitly  use  distinguished 
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boundary  (frequency  response)  information.  However,  the  sequential  manner 
of  the  test  does  use  such  information  in  an  implicit  way.  Also  the  test 
of  Theorem  4 may  be  easier  to  implement  since  the  zero  dependence  of 
A(1 , . . . ,z^,0, . . . ,0)  tends  to  cancel  terms  in  the  original  polynomial 
which  effectively  diminishes  the  complexity  of  the  test. 

At  this  point  we  state  and  prove  the  final  theorem  of  the  section. 

This  theorem  combines  the  "n"  Nyquist  plots  of  Theorem  3 into  a single 
Nyquist  plot. 

Theorem  5:  Let  A be  a polynomial  mapping  £n  into  £n.  Then  A has  no 
zeros  in  Pn  if  and  only  if 

(i)  A has  no  zeros  on  Tn,  and 

(ii)  The  Nyquist  plot  for  the  single  variable  function  A(z,z,...,z) 
does  not  encircle  zero. 

Before  giving  the  formal  proof  of  this  theorem,  we  will  sketch  its  de- 
rivation as  a corollary  to  Theorem  3.  First  observe  that  since  an  analytic 
function  is  orientation  preserving,  [15],  the  degree  of  each  of  the  single 
variable  Nyquist  plots,  Aor(l ,1 , . . . ,1 , • ,1 , . . . ,1 ) of  Theorem  3,  is  non- 
negative. Moreover,  the  degree  function  associated  with  the  set  of 
closed  curves  is  an  isomorphism  from  ir(C  - {0})  onto  the  additive  group 
of  integers.  Thus  the  degree  of  the  concatenation  (the  group  operation 
in  ir(C  - {0}))  of  closed  curves,  will  be  the  sum  of  the  degrees  of  the 
individual  curves.  Furthermore,  the  sum  of  non-negative  integers  is 
zero  if  and  only  if  each  integer  is  zero.  This  implies  that  the 
second  hypothesis  of  Theorem  3 holds  if  and  only  if  the  single  Nyquist 
plot,  obtained  by  concatenating  the  "n"  Nyquist  plots  of  Theorem  3, 
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does  not  encircle  zero.  Consequently  Theorem  3 reduces  to  verifying 
that  A has  no  zeros  on  Tn  and  checking  their  encirclement  of  zero  by 
the  single  Nyquist  plot: 

CAorC- .1 1 )3*CAorCl ,l l )1*  • ••  *[Aor(i 1 .-)] 

= Ao[r(-,l l)Ml,«,l,...»D*  ...  *rU... .,!,•)]  (13) 

where  the  equality  of  equation  [13]  is  due  to  the  fact  [11]  that  com- 
position distributes  over  concatenation. 

Although  equation  13  reduces  the  multivariable  stability  test  to 
a single  Nyquist  plot,  this  plot  is  just  the  concatenation  of  the  n 
plots  of  Theorem  3.  The  resultant  test,  then,  is  no  easier  to  implement 
than  the  original  test.  Fortunately,  this  Nyquist  plot  is  homotopic 
in  C - {0}  to  the  Nyquist  plot  for  the  single  variable  function 
A(z,z,. . . ,z),  obtained  by  setting  each  of  the  dependent  variables 
of  A equal  to  one  another.  To  verify  this  contension,  first  observe 
that  this  single  variable  Nyquist  plot  is  equal  to  the  Nyquist  plot 

Aor(* where  r(*,* *):T^ — Cn  by  taking  the  point  o in  T* 

to  (ct, a,..., a)  in  Cn.  Now,  A has  no  zeros  on  Tn.  Hence  a maps  Tn 
continuously  to  C - {0};  The  Nyquist  plot  Aor(.*, *,...,•)  will 
thus  be  homotopic  in  C - {0}  to  the  Nyquist  plot  of  equation  13 
provided  their  corresponding  Nyquist  contours  are  homotopic  in  Tn. 

This  is,  indeed,  the  case.  However,  the  required  homotopy  is  extremely 
complex.  As  such,  rather  than  wading  through  the  details,  we  will 
simply  sketch  the  required  homotopy  in  the  two  variable  case.  Then 
we  proceed  to  an  alternate  proof  of  the  Theorem  based  on  a known, 
but  non-intuitive,  theorem  of  functions  of  several  complex  variables. 


To  sketch  the  required  homotopy  in  the  two  variable  case,  represent 
2 

the  torus,  T , as  a square  on  the  plane.  Topologically  identify 

opposite  sides  of  the  squares.  Figure  2 illustrate  such  a square. 

The  point  (z-j ,z2) , |z1 1 * |z2|  s 1,  in  T , corresponds  to  the  point 

(e-j,e2)  on  the  plane  where  6-|  = argCzj)  and  @2=  arg(z2)  i.e.  (z-j ,z2)  = 

(e10l,e102).  In  otherwords  the  upper  and  lower  boundaries  of  the  square 

• • 

represent  the  same  line  on  the  torus  since  elir  = e~llr  and  similarly  for 
the  right  and  left  boundaries  of  the  square.  Moreover,  all  four 

corners  of  the  square  represent  the  same  point,  (.-1,-1).  In  the 

sketch  of  Figure  2,  the  Nyquist  contour  r(* ,1 )*r(l »•)  of  equation  [13] 

corresponds  to  the  curve,  number  1,  which  starts  at  (e-|,e2)  = (0,0)  in 

the  center  of  the  square  traveling  vertically  to  the  top  of  the  square. 

It  then  goes  from  the  bottom  of  the  square  vertically  back  to  the 

center,  from  the  center  of  the  square  curve  1 then  passes  horizontally 

to  the  right  hand  boundary,  and  finally  it  returns  from  the  left  hand 

boundary  of  the  square  back  to  the  center.  Since  the  upper  and  lower 

boundaries  of  the  square  are  identified,  when  the  curve  "jumps"  from 

the  upper  to  lower  boundary,  it  remains  continuous  (think  of  the  square 

being  rolled  up  into  a cylinder  with  the  upper  and  lower  boundaries 

glued  together).  Similarly  for  the  "jump"  from  the  right  to  left 

boundary.  Of  course,  curve  1 is  closed  since  it  starts  and  ends  at  the 

same  point. 

The  Nyquist  contour  r(*,*)  is  represented  in  Figure  2 by  curve  5, 
which  starts  at  the  center  of  the  square,  goes  diagionally  to  the  upper 
right  hand  corner  of  the  square  and  then  "jumps"  to  the  lower  left  hand 
corner  of  the  square  from  which  point  it  returns  diagionally  back  to 
the  center.  As  before  the  curve  is  continuous  and  closed. 
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The  required  homotopy  between  curves  1 and  5 is  indicated  on  Figure  2 
by  the  three  intermediary  curves  numbered  2,  3,  and  4.  As  before  these 
curves  are  continuous  since  the  upper  and  lower  boundaries  and  the  left 
and  right  boundaries  of  the  square  are  identified.  Also  all  the  inter- 
mediary curves  begin  and  end  at  the  base  point  (q-|  » ©2 ) * (0,0).  The 
continuity  of  the  intermediary  curves  is  illustrated  graphically  in 
Figure  3 wherein  we  have  redrawn  Figure  2 with  the  point  (e^^)  = 

(+TT,  jj)  taken  as  the  center  point.  In  this  representation  it  is  clear 
that  curves  2,  3,  and  4 are  continuous  and  converge  to  curve  5. 

• 

Although  the  homotopy  required  to  complete  our  proof  is  neatly 
illustrated  in  Figure  2,  its  explicit  mathematical  discretion  is  by 
no  means  simple,  even  for  the  two  variable  case.  Consequent  rather  than 
formalizing  the  tedious  details  of  the  required  n-variable  homotopy  we 
construct  an  alternate  proof  of  Theorem  5 based  on  a theorem  of  several 
complex  variables.  Since  the  theorem  is  applicable  to  analytic  functions 
as  well  as  polynomial:  this  proof  will  also  allow  us  to  extend  theorem 
5 to  the  case  of  meramorphic  transfer  functions  in  several  complex 


variables. 


Lemma  5:  Let  f = (f^ ,f2»* • • ,fn)  be  a continuous  function  mapping  P^ 

to  Cn,n  >_  2,  such  that  f(T^)e;Tn  and  each  of  its  coordinate  functions, 

f.- , have  positive  degree  when  viewed  as  functions  from  T1  to  C,  f,-  1:T1-*C. 

1 t1 

Then  for  any  analytic  function  g:Cn-»C,  g has  a zero  in  Pn  if  and  only 
if  g has  a zero  in  TnUf(P^).  The  theorem  appears  on  page  87  of  reference 


17  and  its  proof  will  not  be  repeated  here.  In  essence  the  theorem  yields 
an  entire  family  of  n-dimensional  Hurwitz-like  tests  (since  Tn  is  n-dimen- 
sional  and  f(P^)  is  2-dimensional)  one  for  each  f satisfying  the  hypotheses 
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of  the  theorem. 

Proof  of  Theorem  5:  To  prove  theorem  5,  we  apply  Lemma  5 with  f defined 

by  f(z)  = (z,z z).  Since  each  coordinate  function  f.(z)  = z is  the 

identity  map  the  coordinates  all  have  degree  one  the  hypotheses  of 
Theorem  5 are  satisfied.  As  such,  the  polynomial  A:Cn  -*■  C (an  arbitrary 
analytic  function,  g:Cn  C,  could  be  used  with  equal  validity)  has  a 
zero  in  Pn  if  and  only  if  it  has  a zero  in  TnUf(p^).  Now, 


f(p') 


{(z,z,.. .,z)  in  Cn 


|z|<l} 


04) 


is  just  a polydisk  in  one  variable  embedded  in  Cn.  As  such,  via  the 
classical  single  variable  Nyquist  criterion  the  existence  of  zeros  of 
A in  f(P^)  may  be  determined  using  the  Nyquist  plot  Aor(* ,•)  whose 
Nyquist  contour  r(*,*,...,*)  follows  the  boundary  of  fCP1).  This  Nyquist 
plot  is,  however,  just  the  classical  Nyquist  plot  for  the  single  variable 
function  A(z,z, . . . ,z) . Thus,  if  we  check  to  see  if  A has  no  zeros  on 
Tn  and  that  the  Nyquist  plot  for  A(z,z,...,z)  does  not  encircle  zero 
we  are  assured  that  A has  no  zeros  in  TnUf(P^)  and  thus  by  Theorem  6 
we  are  assured  that  A has  no  zeros  in  Pn  as  was  to  be  shown. 


Examples 

Example  1 

Consider  the  six-variable  fourth  order  polynomial 


A(2l,z2,z3,z4,z5,z6)  = 10ziz2z3z4zg  + 2z^z2  + z^\z\  + zj?  + 3 (15) 


For  which  the  image  of  A restricted  to  Tn  is  plotted  in  Figure  4.  Since 
zero  is  not  in  the  image  condition  A.  of  the  theorem  is  satisfied  and 
we  may  proceed  to  check  condition  B.  This  requires  that  we  test  the 


Nyquist  plots  for  the  six  one  variable  functions 


A(z-j  ,1 ,1 ,1 ,1 ,1 ) = llz-j  + 2z^  + 4 Q6.) 

A(1  ,Z2»1  »1 ,1 ,1 ) = 12z£  5 (_1 7 • ) 

A(1 ,1  »Z3»1 ,1 ,1 ) = lOz^  z3  ^6  (lg. ) 

A(1 ,1 ,1 ,z^,l ,1 ) - 10z^  + z^  + 6 (lg .) 

A(l,l,l,l,z5,l)  » 10z§  + 7 (20.) 

and 

A(l,l  ,1,1,1  ,z6)  » Zg+  16  (21.) 


for  encirclements  of  zero.  The  resultant  plots  are  sketched  in  Figures 
5a  through  5f.  where  we  see  that  five  of  the  six  plots  encircle  zero. 
As  such,  the  system  is  unstable. 
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Example  2:  In  this  case  we  apply  Theorem  4 to  the  same  polynomial  of 
Example  1.  Again  the  image  of  Tn  under  A does  not  contain  zero  so  we 
must  perform  the  sequential  tests  outlined  in  Theorem  4. 

3 

Step  1 : Consider  the  Nyquist  plot  of  A(1 . ,1 ,Zg)  = Zg  + 16 


Figure  C.  Nyquist  plot  of  A(l,...,Zg)  does  not  encircle  zero. 

Step  2:  Since  the  first  Nyquist  plot  fails  to  encircle  zero,  consider 

2 

the  Nyquist  plot  of  A(1 ,1 ,Zg,0)  = 10zg  + 6.  Clearly  this  Nyquist 
plot  encircles  so  the  filter  is  unstable. 


Figure  7.  Nyquist  plot  of  A( 1 .... ,1 ,Zg ,0) . 
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Example  3:  In  this  we  apply  Theorem  5 to  the  same  polynomial  of  the 
previous  examples.  Notice  that  in  this  case  the  curve  is  of  higher 
order  and  appears  complicated.  There  also  may  be  some  numerical  problems  in 
obtaining  an  accurate  Nyquist  plot  of  this  curve. 

Clearly  A(Tn)  f 0 as  in  the  previous  examples.  Thus  consider  the 


Nyquist  plot  of  AUp...^)  = z7  + 10z6  + 3z3  + 3 where  we  have  set 


Observe  that  the  number  of  times  the  curve  of  Figure  8 encircles 
zero  equals  the  number  of  times  all  the  curves  of  Figure  5,  taken  to- 
gether, encircle  zero. 

Lastly  it  is  interesting  to  wonder  at  the  usefulness  of  these 


plots  for  the  design  engineer.  The  authors  believe  that  in  time  these 
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plots  will  be  shown  to  supply  a large  amount  of  information  on  the  behavior 
of  a system. 

Acknowledgement:  The  authors  wish  to  acknowledge  Rodney  Trotter  who 
ran  numerous  computer  programs  verifying  the  Nyquist  tests  herein. 
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THREE  GRAPHICAL  TESTS  FOR  THE  STABILITY  OF 
MULTIDIMENSIONAL  DIGITAL  FILTERS5 

by 

R.  DeCarlo,  R.  Saeks,  and  J.  Murray 
Abstract 

This  paper  discusses  three  graphical  tests  for  determining  the 
stability  of  multidimensional  digital  filters  characterized  by  an 
appropriate  transfer  function  in  several  complex  variables.  Each 
test  is  carried  out  as  a finite  number  of  "Nyquist"  plots  in  the  com- 
plex plane. 

a 

Introduction 

Recently  two  of  the  authors  constructed  an  algebraic  topological 
proof  of  the  Nyquist  Criterion  (2)  (3).  The  value  of  this  rather 
sophisticated  approach  has  been  harvested  in  generalizations  to  systems 
characterized  by  transfer  functions  in  several  complex  variables  (1) 

(2)  (3),  in  particular  multidimensional  digital  filters.  Specifically 
the  paper  illustrates  three  graphical  tests,  similar  to  the  classical 
Nyquist  test,  carried  out  in  the  complex  plane,  which  determine  the 
stability  of  a multidimensional  digital  filter  with  a transfer  function 
H(z,,...,zn)  = B(z-|,...,zn)/A(z^,...,zn)  where  z..  are  complex  variables 
and  A and  B are  relatively  prime  polynomials.  The  purpose  of  the  paper 
is  to  consider  these  three  tests  as  applied  to  two  different  examples. 

Background  and  Main  Theorems 

Basic  to  the  theory  is  the  2n  (real)  dimensional  polydisc  (8) 
which  is  the  £n  analog  of  the  unit  disc  of  fc.  Mathematically  the  poly- 
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disc  Pn  is 

Pn  = { (21#...,zn)  in  £n  I zi  | <.  1 , 1-1 n} 

There  are  four  separate  notions  of  boundary  of  the  polydisc  ( 1 ) (8). 

First  is  the  usual  topological  boundary 

Bn  = { (z1 zn)  in  ln  | |z.|  <1,1-  1 . ,n 

and  | z ^ | = 1 for  some  k } 

Second  is  the  distinguished  boundary 

Tn  =■  { (z1 zn)  in  en  |z.|  - 1,  i = 1 n} 

Tn  serves  as  the  multidimensional  analog  of  the  ju-axis.  In  particular 
the  frequency  response  (7)  of  a digital  filter  is  the  evaluation  of  its 
transfer  function  over  Tn. 

Thirdly  we  have 

Mn  * { (zj,...,zn)  in  0n  |z.|  = 1,  i = 1 , k-1 , k*l , . . . ,n; 

izki  i' 1 

where  k ranges  from  1 through  n.  This  is  a boundary  notion  in  the  sense 
that  n-1  coordinates  take  on  extremal  values.  Finally,  the  last  notion 
of  "boundary"  is 

Hn  = { (Zl zn)  is  en|  |z.|  =1,  i * 1 , . . . , k-1 ; 

|zkl  < li  z,  =0,  i = k + 1 . ,n> 

where  again  k varies  from  1 to  n.  The  importance  of  this  concept  was 
first  noted  by  Huang  (5).  Later  it  was  generalized  in  (7). 
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With  these  notions  of  boundary  one  may  prove  the  following  Theorem. 
The  proof,  however,  of  the  following  equivalences  is  found  in  the  ref- 
erences (1)  (2)  (5)  (7)  (9). 

Theorem  1:  Let  a causal  multidimensional  digital  filter  be  character- 
ized by  a rational  transfer  function  in  several  complex  variables.  Assume 
the  numerator  and  denominator  polynomials  are  relatively  prime.  Then  the 
following  are  equivalent  stability  conditions: 

(i)  the  pole  set  of  the  transfer  function  has  a null  inter- 
section with  Pn 

(ii)  the  pole  set  of  the  transfer  function  has  a null  inter- 

« 

section  with  Bn 

( i i 1 ) The  pole  set  of  the  transfer  function  has  a null  inter- 
section with  Hn 

(iv)  the  pole  set  of  the  transfer  function  has  a null  inter- 
section with  Mn. 

The  trouble  with  these  conditions  is  that  the  actual  test  is  carried 
out  in  a higher  dimensional  space.  For  example  Pn  is  2n  dimensional,  Bn 
is  (2n-l)  dimensional,  while  Hn  and  Mn  are  both  (n+1)  dimensional.  In- 
tuitively, the  equivalences  of  this  theorem  follow  because  the  pole  set 
of  a rational  function  in  several  complex  variables  is  an  infinite  con- 
tinuum which  must  intersect  the  different  boundaries  of  the  polydisc  if 
it  intersects  the  polydisc  at  all. 

With  the  intuition  gained  in  (2)  (3)  (7)  the  authors  were  able  to 
simplify  these  results  to  graphical  tests  in  the  complex  plane.  The 
following  three  theorems  are  the  fruit  of  this  endeavor.  Before  stating 
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these  theorems,  one  final  definition  is  in  order. 

* The  "Nyquist  plot"  of  a polynomial,  f('),  in  one  complex  variable 

is  defined  to  be  the  image  of  under  the  map  f(')  where  is  the  unit 
circle  of  the  complex  plane.  With  this  in  hand,  we  have  the  following 
three  tests.  Again  the  proofs  can  be  found  in  the  references  (1)  (2) 

(3)  (9). 

Theorem  2:  Let  A be  the  denominator  polynomial  of  a multidimensional 
digital  filter  as  characterized  in  Theorem  1.  A is  a polynomial  mapping 
£n  to  l.  Then  A has  no  zeros  in  Pn  (i.e.  the  filter  is  stable)  if  and 
only  if 

(i)  A has  no  zeros  on  Tn,  and 

(ii)  The  Nyquist  plots  for  the  single  variable  functions 

A(1 1 .z^.l . . ,1 ) k * l,...,n  do  not  encircle  zero. 

Theroem  3:  Let  A be  as  above.  Then  A has  no  zeros  in  Pn  (i.e.  the  filter 
is  stable)  if  and  only  if 

(i)  A has  no  zeros  on  Tn,  and 

(ii)  The  Nyquist  plots  for  the  single  variable  functions 

A(1 ..,1  ,z^,0,...,0)  k = l,...,n,  do  not  encircle  zero. 
Theroem  4:  Let  A be  as  above.  Then  A has  no  zeros  in  Pn  (i.e.  the  filter 
is  stable)  if  and  only  if 

(i)  A has  no  zeros  on  Tn,  and 

(ii)  The  Nyquist  plot  for  the  single  variable  function  A(z,z,...,z) 
does  not  encircle  zero 

Each  of  these  tests  has  essentially  the  same  two  parts.  First  one  per- 
forms the  appropriate  encirclement  test(s);  if  zero  is  not  encircled,  one 
then  proceeds  to  check  the  image  of  the  distinguished  boundary. 
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This  order  of  testing  (encirclement  first,  then  frequency  response)  seems 
in  most  cases  to  be  preferable  to  the  reverse  order,  since  much  less  com- 
putation is  involved  in  the  encirclement  tests;  however,  in  cases  where 
the  frequency  response  is  known  a priori,  or  must  be  plotted  in  any  case, 
the  order  is  immaterial. 

It  might  appear  that  the  third  test  (Theorem  4)  is  the  best,  since 
it  involves  only  one  encirclement  test;  however,  in  many  cases  the  rel- 
ative complexity  of  the  polynomial  A(z,z,...,z)  will  more  than  offset 
this  advantage.  Similarly,  in  many  cases.  Theorem  3 may  be  much  easier 
to  apply  than  Theorem  2.  (This  is  illustrated  in  the  first  example) 
Theorems  2 and  4,  however, do  have  two  advantages.  The  first  is  mainly 
philosophical;  these  Theorems  give  a test  for  stability  purely  in  terms 
of  the  frequency  response  of  the  function  A,  which  corresponds  closely 
with  the  idea  of  the  Nyquist  criterion  in  one  variable.  The  second 
advantage  is  that  by  filling  in  the  interior  of  the  encirclement  plot(s) 
and  taking  this  region  together  with  the  image  of  the  distinguished 
boundary,  one  obtains  the  image  of  the  entire  polydisc,  from  which  one 
can  get  an  accurate  idea  of  stability  margins.  (The  point  here  is  that  we 
have  found  the  image  of  a 2n-dimensional  set--the  polydisc— by  plotting 
an  n-dimensional  set  and  a 1-dimensional  set). 

EXAMPLES 

In  this  section  we  apply  each  of  the  above  tests  to  two  examples. 
Exampl e 1 : 


A(z1,z2)  = 5/4  z*  z|  + 1/2  z^2  + 1/2  zj  z?  + 3 z1  zj;  + 3 z*  z2 
" zi  " z2  + ^ zlz2  ~ ^ Z1  ” ^ z2  + 1 


37 


In  this  case,  we  have 

A(z,l ) = 1/2  z3  + 13/4  z2  + 9/2  z - 2 

A(l,z)  = 1/2  z3  + 13/4  z2  + 9/2  z - 2. 

These  polynomials  are  identical;  the  image  of  the  unit  circle  being  given 
in  Fig.  1(a).* 

Since  this  curve  encircles  the  origin,  we  deduce  immediately  that  the 
filter  is  unstable;  for  purposes  of  illustration,  we  will  carry  out  the 
other  tests. 

A(z,0)  = -z2  - 2z  +1 

A(l,z)  is  as  before  (Fig.  1(a));  A(z  ,0)  is  plotted  (for  z = e19) 
in  Fig.  1(b).  Again,  either  plot  suffices  to  verify  instability,  and 
clearly  A(z,0)  gives  the  simpler  test. 

To  apply  the  third  test,  we  calculate 

A(z,z)  = 9/4  z4  + 6 z3  + z2  - 4z  + 1 

and  the  image  of  the  unit  circle  under  this  mapping  is  plotted  in  Fig.  1(c). 
The  relative  complexity  is  apparent;  however,  it  again  verifies  instability. 
Finally,  we  plot  the  image  of  the  distinguished  boundary  in  Fig.  1(d); 
it  can  be  seen  that  it  does  not  include  the  origin,  although  it  does  in 
some  sense  "encircle"  it. 

The  second  example  shows  that  this  last  kind  of  "encirclement"  is 
irrelevant;  nothing  can  be  deduced  from  it. 


♦(This  illustrates  the  obvious  fact  that  if  the  polynomial  is  symmetric 


in  Zj,...fz  , then  the  n plots  in  Theorem  2 in  fact  reduce  to  1 plot— 
usually  simpler  than  the  plot  in  Theorem  4.  Such  symmetry  is  quite  common). 


Example  2: 

A(zrz2)  = (Zl  + 2)3  (z2  + 2)3 
As  before,  the  plots  for  A(z,l)  and  A(l,z)  are  identical 
A(z,l ) = 27(z  + 2)3 
A(l,z)  = 27(z  + 2)3; 

This  plot  is  given  in  Fig.  2(a);  it  does  not  encircle  0. 

In  this  case,  the  plot  for  A(z,0)  = 8(z  + 2)  differs  from  the  previous 
plot  only  by  a scale  factor;  we  do  not  draw  it  separately. 

Finally,  the  plot  for  A(z,z)  = (z  + 2)5 6  is  given  in  Fig.  2(b); 
again  it  does  not  encircle  the  origin. 

Thus,  in  order  to  determine  stability  in  this  case,  it  is  necessary 
to  plot  the  image  of  the  distinguished  boundary;  this  is  done  in  Fig.  2(c). 
Since  this  image  does  not  contain  the  origin  (although  it  does  surround  it), 
we  conclude  that  the  filter  is  stable.  This  of  course  is  obvious  analytic- 
ally; the  present  example  is  merely  to  illustrate  the  tests. 

Note:  Because  of  the  magnitudes  of  the  numbers  involved.  Figs. 2(a)  - 2(c) 
are  not  drawn  to  scale. 
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A NYQUIST-LIKE  TEST  FOR  THE  STABILITY 
OF  TWO  DIMENSIONAL  DIGITAL  FILTERS* 

by 

R.  DeCarlo,  R.  Saeks,  and  J.  Murray 

ABSTRACT : This  paper  constructs  a Nyquist-like  test  for  the  stability 
o^  two  dimensional  digital  filters.  The  test  takes  the  form  of  a con- 
tinuum of  classical  one  variable  Nyquist  plots  parameterized  by  the 
elements  of  the  unit  circle  of  the  complex  plane.  Since  the  parameter 
space  is  compact  the  test  can  be  accurately  approximated  by  a finite 
number  of  classical  one  variable  Nyquist  plots  and  is  therefore  readily 
implemented  on  a computer. 

A two  dimensional  digital  filter  is  characterized  by  a rational 
transfer  function  in  two  complex  variables 

SUj.Zg) 

1#  A (Zj.Zg) 

where  A (Zj,z2)  and  B(z^,z2)  are  relatively  prime  polynomials  in  z^  and 
z2-  For  the  purpose  of  this  paper  we  say  that  the  digital  filter  is 
stable  if  Afz^Zg)  f 0 for  Iz^  <_  1 and  |z2|  £ 1.  This  structural  stability  con- 
dition implies  that  the  filter  is  bounded-input  bounded  output  stable 
though  as  recently  shown  by  Goodman^  the  condition  is  actually  slightly 
stronger.  Huang  showed  that  this  4-dimensional  stability  condition  was 
actually  equivalent  to  the  3-dimensional  condition  that  A^.Zg)  t 0 for 
|z^|  = 1 and  |z2l  ± 1 or  |z^|  <_  1 and  z2  = 0 which  we  use  as  the  basis 
of  our  theory. 

The  key  to  the  formulation  of  our  Nyquist-like  theory  is  the  obser- 
vation that  from  an  abstract  analytic  function  point  of  view  the  classical 
one  variable  Nyquist  plot  is  simply  a method  of  determining  whether  or  not 
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an  analytic  function  in  one  variable  has  zeros  in  an  appropriate  region 
by  plotting  the  image  of  the  function  on  the  boundary  of  the  region.  To 
obtain  a Nyquist  theory  in  two  variables  we  therefore  decompose  the  region 
of  C in  which  A(z1,z2)  is  forbidden  to  have  zeros  by  Huang's  theorem  in- 
to the  union  of  a family  of  one  variable  regions  to  which  the  classical 

Nyquist  theorem  applies.  More  precisely,  for  each  complex  number  a of 

2 

magnitude  one  we  define  the  disk  Do  in  C by 

2.  0a  = {(e1a  ,z2)  ; |z2|  < 1} 

and  we  define  the  disk  D„  by 

o 

3.  Dq  = {(zr0)  ; |zx|  < 1} 

Now,  Huang's  theorem  may  be  restated  as  "the  digital  filter  is  stable  if 
and  only  if  A(z^,z2)  has  no  zeros  in  the  disks  Dq  and  |aj  * 1". 
Observing  that  each  disk  is  fixed  in  one  of  its  variables  the  polynomial 
A(z1,z2)  restricted  to  any  of  the  above  defined  disk's  is  an  analytic 
function  of  one  variable  and  hence  the  classical  Nyquist  test  can  be  used 
to  check  for  zeros  within  the  disk.  In  particular,  A(z^,z2)  has  zeros  in 
the  disk  if  and  only  if  the  Nyquist  plot  for  the  one  variable  function 
A(el0t.z2)  does  not  equal  or  encircle  zero.  Similarly,  A (z^,z2)  has  no 
zeros  in  the  disk  DQ  if  and  only  if  the  Nyquist  plot  for  the  one  variable 
function  A(z1>0)  does  not  equal  or  encircle  zero.  Combining  these  obser- 
vations we  obtain  the  following  stability  theorem. 

Theorem:  A digital  filter  characterized  by  the  two  variable  transfer 


function 


B(zltz2) 


where  A(z^,z2)  and  B(z^,z2)  are  relatively  prime  polynomials  in  two 
variables  is  stable  (structurally  stable)  if  and  only  if  the  Nyquist 
plots  for  the  family  of  one  variable  functions 

A(ela,z2)  ; |a|  « 1 
and 

AUj.O) 

do  not  equal  or  encircle  zero. 

Although  the  theorem  formally  implies  that  one  check  a continuum 
of  Nyquist  plots  parameterized  by  the  complex  numbers  of  magnitude  one, 
in  fact,  since  this  set  of  numbers  is  compact  one  can  obtain  a test  with 
arbitrarily  good  resolution  using  only  a finite  number  of  plots.  Indeed, 
in  a somewhat  different  context  the  authors  have  shown  that  a similar  con- 
tinuum of  Nyquist  plots  can  actually  be  reduced  to  a single  plot  without 

2 

inducing  any  error  into  the  stability  test.  The  following  examples  are 
based  on  a finite  approximation  to  the  continuum  of  plots  required  by 
the  theorem. 

Example  1:  Let  the  transfer  function  of  a digital  filter  be 


4. 


HCzj.Zg) 


1 

1 + .25ZJ  + .25z2 


b(Zi»z2) 
A(z1,z2)  * 


Step  1:  Draw  the  Nyquist  plot  for  A(z^,0).  This  curve,  shown  in  Figure 
1,  does  not  encircle  zero.  So  we  proceed  to  the  next  step  as  outlined 
in  the  theorem.  iu»/f 


Figure  1:  Nyquist  plot  of  A^.O)  for  example  1. 


Step  2:  Now  consider  the  family  of  Nyquist  plots  for  the  functions 
A(e1a,z^)  ; |a[  = 1.  This  family  of  curves  does  not  encircle  "0" 
as  indicated  in  Figure  2.  Thus  the  filter  is  stable. 


Z/=  i 


z,=.7l  -t-i  .71 


Figure  2:  Nyquist  plots  for  A(e  .Zg)  for  example  2. 


Example  2:  Now  consider  the  filter  whose  transfer  function  is 


t is  illu 


• If# 
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At  this  point  no  decision  can  be  made  so  proceed  to  step  2. 
Step  2:  Consider  the  family  of  functions  A(e1ci,z2)  , |a|  = 1. 
Nyquist  plots  for  some  of  these  functions  are  shown  in  Figure  4. 
They  indicate  that  the  filter  is  indeed  unstable. 


Figure  4:  Nyquist  plots  for  A(e1a,z2)  for  example  2. 
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Another  Proof  and  a Sharpening  of  Huang's  Theorem 

By 

John  Murray 

ABSTRACT : Two  further  proofs  of  Huang's  Theorem  on  the  zeros  of  analytic 
functions  of  two  variables  are  given;  the  first  is  similar  to  previous 
proofs,  but  is  made  shorter  by  the  use  of  a known  maximum  modulus  principle; 
the  second  is  completely  different,  using  a theorem  of  Rudin  which  actually 
gives  a sharper  result  than  Huang's.  Finally  it  is  indicated  how  a corres- 
pondingly sharper  result  may  be  obtained  in  higher  dimensions. 

In  [5]  Huang  stated  and  gave  an  incomplete  proof  of  the  following 
theorem: 

THEOREM: 

A two-variable  polynomial  P(Z^,Z2)  has  a zero  in  the  polydisc  (J2  = 

{ (Zj ,Z2 ) |Zj|  _<  1,  | Zg | 1 1)  if  and  only  if  it  has  a zero  in  the  set 
H2  = {(zrz2)  | |Zj|  = 1,  |Z2|  <1  }U  { (Zj ,0)  | |Z1i  <1  } 

Two  proofs  of  the  Theorem  have  appeared  £2,3],  both  appealing  to  the 
maximum  modulus  principle;  the  first  proof  here  is  similar  but  shorter  since 
.it  uses  a convenient  known  version  of  the  maximum  modulus  principle. 

We  will  use  the  following  notation: 

U = {Z|  |Z|  <1} 

U = {Z|  |Z|  <1} 

T = f z ! |Z!  = 1} 

Un  = Uxllx...xU  (n  times) 

If1  = UxUx...xU  (n  times) 

Tn  = TxTx...xT  (n  times) 


= {(ZrZ2) 
= {(Zr0) 
Then  H?  = H£  UH£. 


IZil  = i. 
IZjl  11} 


2 1 - 


< 1} 
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in  W such  that  Z ^ = 0,  i.e.,  there  is  a zero  in  HjJ 


Proof  I: 

a a __ 

Assume  PfZ^.Zg)  has  a zero  (Z^,  Z2)  in  Ug.  Let  V denote  the  connected 
component  of  the  zero-set  of  which  contains  (Z^.Z^)  and  let  W = 

vn{(Z1,Z2)  |Zj|  <1}.  If  (z^,z|)  i W,  (Z*,z|)  e H^  = H2;  we  can  therefore 
assume  that  (Z^.Z^)  e W f 0.  We  consider  the  function  1/Z2  restricted 
to  the  analytic  set  W;  if  this  function  is  not  analytic,  there  is  a point 

If  the  function  is 

analytic,  then  by  the  maximum  modulus  principle  [4,p.  106]  |1/Z2|  can  not 
have  a strict  maximum  in  W.  Thus,  since  | Z^ | 1 1 and  (Zj,  Z^)  e W,  there 
is  a point  (Z^.Z^)  e W H{ (Z^^ , Zg ) | Z ^ | = 1}  such  that  | Z^ I ± 1.  By  continuity, 
P(zJ,  z|)  = 0;  thus  (Zj ,z|)  is  a zero  of  P in  H£.  Q.E.D. 

Before  proceeding  to  the  second  proof,  one  further  definition  is  needed: 

If  y:  U -*•  U is  continuous,  by  Ind  >PoE  we  will  mean  the  usual  winding  number 

O • A 

of  the  closed  curve  f(e<:iTj,  0 <_  t <_  1,  about  the  origin.  (In  what  follows, 
this  curve  will  never  go  through  the  origin). 

Proof  II: 

We  will  use  the  following  theorem  of  Rudin  [6,p.  87] 

If  i = (ij.Jg,.  ••  ,$n):  U -►  U11  is  a continuous  mapping  such  that  i(T)  =Tn 
and  Ind  5i  oE  > 0 for  all  i,  then  f (Z^.Zg, . . .Zn)  has  a zero  in  U11  if  and  only 
if  f(Zj,...,Zn)  has  a zero  in  the  set  Tn  U 4(U).  (Here  f is  any  function 
analytic  on  Un  and  continuous  on  U?) 

To  apply  this  theorem,  we  take  n = 2,  and  define  4(Z)  as  follows: 


*(Z)  = ( 


(2Z.0) 


|Z|  < >1 


(Z/|Z|,  (2 | Z | -1)Z)  , 1/2  < | Z | < 1. 


4 is  clearly  continuous,  and 

i.(e2irit)  » e2iriJ  i = 1,2, 
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so  4(T)  = T2,  and  Ind  &.  oE  = 1,  i = 1,2. 

Thus  if  has  a zero  in  ir,  it  has  a zero  in  T U 4 (U);  this  set 

is  clearly  a subset  of  Hg,  and  the  result  follows.  Q.E.D. 

We  note  the  following:  Firstly,  we  have  proved  only  one  direction 

of  the  implication:  the  other  is  trivial.  Secondly,  both  of  the  above 

2 —2 

proofs  apply  to  functions  analytic  on  U and  continuous  on  U and  not 
merely  to  polynomials.  Thirdly,  proof  II  yields  a considerably  sharper 
result  than  Proof  I (or  any  previous  proofs  of  Huang's  Theorem)  in  that 
the  set  which  must  be  tested  for  zeros  is  a proper  subset  of  H;  in  fact 
dim  (T  U 4 (H))  = 2,  while  dim  H = 3.  Fourthly,  a correspondingly 
sharper  version  of  the  higher-dimensional  extension  of  Huang's  Theorem 

[1]  can  be  obtained  by  defining 
4(Z):  U-U"  by 

fi(z)  = (-[7p • • • »]7p  ^ (nUI  " k + o....»o),  — <.  | z { <_£ 

t 

kth  entry 

for  1 ± k <_  n. 

Finally,  other  choices  of  4 will  yield  other,  possibly  simpler,  tests; 
this  is  explored  in  [7]. 
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ABSTRACT 

Two  sets  of  necessary  and  sufficient  conditions  for  the  existence 
of  a polynomial  spectral  factorization  of  a given  polynomial  are  proved. 
These  conditions  are  shown  to  be  necessary  and  sufficient  for  the 
existence  of  a rational  spectal  factorization  of  the  spectral  function 
of  any  rational  function  and  necessary  but  not  sufficient  for  the  ex- 
istence of  a spectral  factorization  of  the  rational  function  itself. 

Here  the  term  "spectral  factorization"  is  to  be  understood  as  a factor- 
ization into  the  product  of  a function  without  poles  and  zeros  in  the 
open  unit  polydisc,  and  a function  without  poles  or  zeros  in  a region 
inverse  to  the  open  unit  polydisc.  These  conditions  are  seen  to  give 
extremely  severe  restrictions  on  the  spectral  function  of  the  given 
polynomial  or  rational  function,  and  hence,  on  the  amplitude  response 
of  any  possible  quarter-plane  purely  recursive  (stable)  digital  filter. 
The  implications  of  these  restrictions  for  the  design  of  minimum-phase 
finite-impulse-response  and  stable  infinite-impulse-response  filters 
are  discussed.  In  particular,  it  is  shown  that  the  difficulties 
which  many  researchers  have  encountered  in  stabilizing  two-dimensional 
infinite-impulse-response  filters  are  inherent  in  the  problem  and  can 
not  be  avoided  by  a refinement  of  technique;  any  method  which  attempts 
to  stabilize  a filter  by  finding  a stable  denominator  polynomial  whose 
spectral  function  matches  (or  closely  approximates)  that  of  a given  un- 
stable denominator  polynomial  must  fail  for  a large  class  of  polynomials. 
This  is  because  there  is  no  stable  polynomial  having  the  same  spectral 
function  as  a polynomial  in  this  class  (nor  even  approximating  it  well) ; 
thus,  it  is  precisely  the  attempt  to  match  the  amplitude  response  of 
the  unstable  filter  which  forces  the  new  filter  to  be  unstable  also. 

This  is  in  sharp  contrast  to  the  situation  in  one  dimension,  where  any 
rational  spectral  function  has  a rational  spectral  factorization. 
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Introduction: 

The  subject  of  two-dimensional  digital  filters  has  received  consi- 
derable attention  of  late:  in  particular,  two-dimensional  spectral  fac- 
torization has  been  treated  in  a number  of  papers  - it  is  considered  in 
great  detail  in  reference  [1],  The  major  problem  which  arises  is  that  in 
general  the  spectral  factors  of  a rational  transfer  function  are  not 
rational:  some  further  processing,  such  as  truncation  and  smoothing,  is 
usually  employed  to  yield  approximate  rational  factors.  It  is,  therefore, 
somewhat  surprising  that  the  class  of  rational  functions  for  which  a ra- 
tional spectral  factorization  exists  does  not  seem  to  have  been  investi- 
gated. In  this  paper,  we  give  two  sets  of  conditions  which  must  be  satisfied 
by  such  functions  (theorems  1 and  3);  a converse  is  given  which  may  be 
applied  to  the  numerator  and  denominator  polynomials  separately.  Now,  the 
polynomial  spectral  factors  (when  they  exist)  of  a given  polynomial  are 
minimum-  and  maximum-phase  polynomials;  conversely,  every  such  polynomial 
gives  rise  to  trivial  spectral  factors.  Motivated  by  this,  we  apply  the  re- 
sults of  theorems  1 and  3 to  the  particular  case  of  minimum-phase  polynomials 
(i.e.,  polynomials  without  zeros  in  the  unit  polydisc). 

In  this  context,  the  main  consequences  of  the  results  of  this  paper  may 
be  broadly  outlined  as  follows: 

i)  A given  polynomial  has  exactly  the  same  amplitude  response  as  a 
minimum-phase  polynomial  if  and  only  if  the  classical  one-variable 
method  (of  factoring  the  original  polynomial  into  a product  of  two 
polynomials  devoid  of  zeros  in  certain  regions)  can  be  applied. 

(This  result  is  in  fact  implicit  in  [1],  but  does  not  appear  to 
have  been  explicitly  stated  in  the  literature).  The  corresponding 
statement  for  minimum-phase,  stable  rational  functions  is  false, 
however. 
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11)  If  the  conditions  given  in  theorems  1 and/or  3 are  not  satisfied, 
then  not  alone  is  there  no  minimum-phase  stable  rational  function 
having  exactly  the  same  amplitude  response  as  the  original;  the 
original  amplitude  response  can  not  even  be  approximated  arbi- 
trarily well  by  minimum-phase  stable  rational  functions.  This 
follows  from  the  fact  that  the  conditions  in  theorems  1 and  3 are 
conditions  on  the  amplitude  response  which  are  preserved  under 
any  reasonatie  kind  of  convergence. 

iii)  The  conditions  in  theorem  3 are  easily  visualized  and  surprisingly 
stringent;  they  require  essentially  that  the  gain  of  the  filter, 
averaged  over  certain  directions  in  the  frequency  plane,  have  no 
variation  in  a perpendicular  direction,  (see  the  discussion  follow- 
ing theorem  3).  This  gives  extremely  severe  restrictions  on  the 
amplitude  response  of  minimum-phase  FIR  filters,  minimum-phase 
stable  IIR  filters,  and  the  denominator  polynomial  of  arbitrary 
stable  IIR  filters. 

iv)  It  has  been  pointed  out  by  Bose  [9]  and  Woods  [10],  and  again  is 
implicit  in  [1],  that  there  exist  purely  recursive  filters  whose  am- 
litude  responses  are  not  realizable  as  the  amplitude  response  of 
any  stable  purely  recursive  filter,  and  that  consequently  any 
stabilization  method  which  attempts  to  match  the  amplitude  re- 
sponse of  the  original  filter  is  doomed  to  failure.  The  restrict- 
ions referred  to  in  iii), above,  reinforce  this  conclusion  and 
identify  the  precise  properties  of  the  examples  in  [9]  and  [10] 
which  make  stabilization  Impossible. 

Definitions  and  Notation: 

Our  notation  will  follow  that  in  [2];  we  repeat  it  here  for  convenience. 

For  simplicity  we  restrict  ourselves  throughout  to  two  dimensions,  although 
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there  does  not  appear  to  be  any  difficulty  in  extending  the  results  to 
higher  dimensions.  Thus  all  functions  are  assumed  throughout  to  be  rational 
functions  of  two  complex  variables  unless  otherwise  stated;  we  further  ex- 
clude the  zero  function.  Two-dimensional  complex  space  will  be  denoted 

2 2 

by  If  ,i.e,  $ s {(Zj,  Z2)|  1\  and  Z2  are  complex  numbers}.  The  open  unit 

2 

polydisc  will  be  denoted  by  U ,i.e., 

U f ( (Zi ,Z2)e  If  | |Zj  | <1  and  |Z2|  <1} 

2 

and  its  closure  will  be  denoted  by  U : 

lf=  {(lu  Z2)e  <f2|  |ZL | <1  and  |Z2|  _<1} 

2 

The  distinguished  boundary  of  the  unit  polydisc  will  be  denoted  by  T : 

T2=  ( ( Z ! , Z2)e  <f2|  | Z ! | = 1 and  |Z2|  = 1} 

The  frequency  response  of  the  filter  whose  transer  function  is  f(ZlfZ2)  is 

2 

simply  the  restriction  of  f to  T . We  will  find  it  convenient  to  denote  this 
restriction  by  f*. 

The  one-dimensional  sets  corresponding  to  the  above  are: 

U = (Z  £ (f  I |Z|  <1} 

U = (Z  £ If  I |Z|  <1} 

T - (Z  £ (I  | Z | =1} 

.2 

We  need  one  further  subset  of  If  : 

V2=  { ( Z ! , Z2 ) e If2  | | Z ! | >1  and  |Z2|  >1}  . 

2 

By  the  Fourier  coefficients  of  a function  h (0i,02)  defined  on  T we 


Li 


mean  the  numbers 
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amn  — 2 


f 2'  h(0„92)  e-J(m9‘+  "9*>d9,  de2. 


Finally,  let  us  state  precisely  what  we  mean  by  the  term  spectral 
factorization.  Several  different  forms  of  spectral  factorization  are  treated 
in  [1];  here  we  will  be  concerned  only  with  the  simplest  form:  if  f is  a 
rational  function,  it  will  be  said  to  have  a (rational,  quarter-plane) 

spectral  factorization  if  f = f\f2  where  fx  and  f2  are  rational  functions, 

2 2 

fi  has  no  poles  or  zeros  in  U , and  f2  has  no  poles  or  zeros  in  V . Several 
conments  are  in  order  concerning  this  definition: 

i)  By  "rational"  we  mean  only  "finite-order";  i.e.,  the  functions 

are  assumed  to  be  expressible  as  the  quotient  of  two  (finite-order) 
polynomials. 

ii)  The  quarter-plane  property  enters  only  in  connection  with  the 

regions  in  which  the  factors  are  assumed  to  be  zero-  and  pole-free; 

2 

in  particular,  if  f has  no  poles  or  indeterminacies  on  T , and  has 
a quarter-plane  spectral  factorization,  then  there  is  a quarter- 
plane  causal,  stable  filter  whose  amplitude  response  is  equal  to 


iii)  It  would  possibly  be  more  natural  to  work  with  IP  and  V2  rather  than 
U2  and  V2  (especially  when  considering  stability).  However,  to  do 
so  would  complicate  the  statements  of  the  theorems  considerably,  and 
it  is  usually  clear  whether  or  not  the  results  will  hold  with  U2  and 

V2  in  place  of  U2  and  V2.  (One  needs  only  to  check  for  zeros  and 

2X 

poles  on  T ).  In  general,  if  the  "closed"  version  is  not  obvious, 
it  is  not  true;  1 - Zx  Z2  will  serve  as  a counterexample  in  all 
such  cases. 
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iv)  To  simplify  the  statements  of  the  theorems,  the  definition  has 

been  given  in  terms  of  the  rational  function  f itself,  rather  than 

2 

the  spectral  function  Jf*|  ; however,  the  conditions  given  in  the 

I .2 

theorems  actually  involve  only  |f*|  . 

, 2 2 

v)  We  note  that  V is  defined  to  be  a subset  of  $ ; thus  the  behaviour 

of  functions  at  infinity  is  irrelevant  to  our  purposes. 

Spectral  Factorization: 

Our  first  criterion  for  the  existence  of  rational  spectral  factors  is 
very  much  in  the  spirit  in  which  spectral  factorization  is  treated  in  [1]; 
it  is  a trivial  consequence  of  theorem  5.4.7  in  [2]. 

Theorem  1 : 

2 

If  a rational  function  f on  f has  a rational  spectral  factorization 
then  the  Fourier  coefficents  a „ of  log  |f*|  are  zero  for  all  pairs  of  in- 
teg ers  (m,n)  such  that  m t o,  n f o,  and  m and  n have  different  signs  - 
that  is,  for  all  integer  points  in  the  second  and  fourth  quadrants  The 
converse  is  true  for  polynomial  f. 

As  mentioned  above,  this  criterion  involves  only  the  absolute  value 
of  f;  it  follows  that  the  existence  of  spectral  factors  imposes  restrictions 
on  the  amplitude  response  of  a two-dimensional  filter  - in  contrast  with 
the  situation  in  one  dimension.  The  above  criterion,  however,  does  not 
present  these  restrictions  in  an  easily  visualized  form  - for  instance,  it 
is  difficult  to  guage  exactly  how  severe  the  restrictions  are.  For  this 
reason,  we  next  present  conditions  which  are  stated  in  terms  of  the  log- 
amplitude  response  itself,  rather  than  its  Fourier  coefficients.  This 
result  takes  an  approach  which  seems  to  differ  substantially  from  those 
previously  known;  it  gives  easily  visualized  necessary  conditions  on  those 
rational  functions  which  admit  a rational  spectral  factorization.  Before 
we  state  this  theorem,  however,  we  first  present  a simple  result  which  will 
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be  used  in  the  proof,  and  is  also  of  separate  interest;  one  of  its  conse- 
quences is  that  when  rational  spectral  factors  exist,  the  usual  one-dimensional 
stabilization  method  (for  unstable  denominator  polynomials)  can  be  used. 

Theorem  2: 

If  the  rational  function  f admits  a rational  spectral  factorization, 
then  there  is  a rational  function  f (with  deg  f <_  deg  f)  such  that 

|f*|  = |f*L 

^ 2 

and  f has  no  poles  or  zeros  in  U . 

Again,  the  converse  holds  for  polynomial  f. 

Thus,  if  the  denominator  polynomial  of  an  unstable  filter  has  polynomial 

spectral  factors,  there  is  a stable  filter  of  at  most  the  same  order  with 

, 2 

the  same  amplitude  response  (provided  the  polynomial  has  no  zeros  on  T ). 

Again  most  of  the  proof  is  contained  in  [2];  we  fill  in  the  details 

p 

here:  suppose  f has  rational  speptral  factors,  then  f ® fj  jj-  where  fi  has  no 

2 2 

poles  or  zeros  in  U and  P and  Q are  polynomials  without  zeros  in  V . 

^ m n 

Let  P = Zi  Z2  P (1/Zj ,1/Z2) , Z2  f 0,  Za  t 0 

where  m is  the  degree  of  P in  Zlf  n is  the  degree  of  P in  Z2,  and  P is  the 
polynomial  whose  coefficients  are  the  complex  conjugates  of  the  coefficients 

of  P.  Clearly  Pisa  polynomial  of  degree  less  than  or  equal  to  the  degree 

'1/ 

of  P,  and  so  is  also  defined  for  Z^O  and  Z2=0.  Now  if  P ( Zx , Z2 ) * 0 for 
Zj  f 0 and  Z2  / 0,  then  F(l/Zlf  1/Z2)  * 0;  this  implies  that  either 

2 

11/ZjI  1 1 or  1 1 /Z2 1 <_  1 (since  P has  no  zeros  in  V ) 

and  so  either 

2 

|Zi|  > 1 or  |Z2|  > 1,  i.e.,  (Zlt  Z2)  * U 


L 


t 


t 
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Thus  the  only  possible  zeros  of  P in  U are  for  lx  = 0 or  Z2  * 0.  But 
by  standard  results  in  the  theory  of  several  complex  variables  [8],  if 
the  zero-set  were  nonempty,  this  would  imply  that  either  lx  or  Z2  was  a 

factor  of  P,  which  is  impossible  by  our  choice  of  m and  n.  Thus  P has 

2 2 

no  zeros  in  U . Finally,  on  T 

'v  m m _ 

|P(ZlaZ2)|  - I Zi  Z2  P (I/Z1.I/Z2) | - |P(Zi.Z2) | = |P(Zx,Z2)|. 
Q is  defined  similarly  and  has  similar  properties.  Then 


clearly  has  the  required  properties. 

Conversely,  suppose  f is  any  polynomial  for  which  there  is  a 

^ 2 

rational  function  f without  poles  or  zeros  in  U such  that 


|f*l  = If*l 

% 

then  f/f  is  rational  and  analytic  in  U , and 


i(f/f)*i  = i 


Thus  by  theorems  5.2.5  and  5.2.6  in  [2],  f/f  = P/Q  where  P and  Q are 

, 2 2 

polynomials,  P has  no  zeros  in  V , and  Q has  no  zeros  in  U . Then 


f = Pf/Q 

gives  a rational  (in  fact,  polynomial)  spectral  factorization  of  f. 


The  Second  Criterion: 

Our  second  set  of  conditions  for  the  existence  of  a rational 


spectral  factorization  is  given  in  the  following: 


A 
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Theorem  3: 

12 

If  a rational  function  f on  (f  admits  a rational  spectral  factori- 
zation, then 

1 

U2;  log  ->)|de 

is  a constant  independent  of  4* , (0  <_  m*  < 2 it) , for  all  integers  m > o and 

i 

n > o. 

Again,  these  conditions  depend  only  on  the  amplitude  response  of  f. 

The  simplest  condition  is  that  for  m = 1 and  n = 1;  it  can  be  easily  visua- 
lized by  drawing  two  adjacent  squares  in  the  0i  02  “ plane  on  which  the 
amplitude  response  is  defined  (the  frequency  response  extends  to  the  entire 
0^2  - plane  by  periodicity),  and  drawing  lines  with  slope  1 and  length 
2-nS?  on  these  squares;  see  figure  1. 

Then  the  condition  for  m = 1,  n = 1 can  be  restated  as:  the  "average" 
amplitude  of  the  function  f along  the  line  is  a constant  - that  is,  it 
is  as  independent  of  the  particular  line  L..  chosen.  ("Average"  here  is  to 
be  understood  as  the  geometric  mean  of  the  amplitude,  or  the  arithmetic 
mean  of  the  log-amplitude).  Alternatively,  we  may  say  that  the  average 
level  of  the  amplitude  over  any  line  of  slope  1 and  of  length  2 is  in- 
dependent of  the  position  of  the  line  in  the  o1e2-plane.  (For  example,  we 

A 

could  vary  the  over  the  dotted  square  in  the  direction  n).  The  condi- 
tions for  higher  m and  n have  a similar  interpretation,  with  a slope  of 
n/m  instead  of  1,  and  length  2ir/m2  + n2  instead  of  2tt/7;  clearly,  if  m and 
n are  not  relatively  prime,  the  corresponding  condition  is  superfluous. 

This  theorem  then  gives  a striking  limitation  on  the  amplitude  re- 
sponse of  a rational  function  which  admits  a rational  spectral  factori- 
zation; even  the  simplest  of  the  conditions  (that  for  n=m*l)  implies  that 

1 J 
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such  a function  can  not  accurately  approximate  an  amplitude  which  has 

A 

large  variations  in  overall  level  in  the  direction  ^ shown  in  figure  1. 
Proof  of  Theorem  3: 

In  view  of  theorem  2,  it  suffices  to  prove  this  under  the  assumption 

2 

that  f has  no  poles  or  zeros  in  U . This  assumption  implies  that  f has  a 

2 

holomorphic  logarithm  in  U . Then,  for  any  integers  m > o,  n > o and  any 
real  number  f, 

log  f(Zm,  Zn  ejv) 

is  a holomorphic  function  of  one  complex  variable  for  Z e U.  Thus 
Re(log  f(Zm,  Zn  ejy)) 

is  a harmonic  function  in  U,  and  so  by  the  mean-value  property  of  harmonic 
functions 

\ Re  (log  f(Zm,  Zn  e3’1'))  de  = Re  (log  f(om,  on  e'3'1')) 

2tt  t 

i.e.,  1 Re  (log  f(e3nie,  e^n9+'t’3))  de  = Re(log  f(o,o)) 

2 it  Jo 

But  Re  log  w = log  |wj  for  w f o,  and  so 

^ ^ log  |f(e3m9,  e'3(n9+4'3)|  de  = log  |f(o,o)| 

and  the  right-hand  side  is  independent  of  f (and,  incidentally,  of  m and 
n also). 

An  obvious  question  which  arises  is  the  extent  to  which  the  converses 
of  these  results  hold.  In  fact,  the  converse  of  Theorem  3 holds  for  poly- 
nomials, and  modified  converses  of  both  Theorems  1 and  3 hold  even  for 
rational  functions.  The  modification  takes  the  following  form:  If  the 
Fourier  coefficients  of  log  |f*|  (where  f is  a rational  function)  vanish 
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for  mn  < o,  then  there  is  a rational  function  f with  rational  spectral 
factors,  (equivalently,  a rational  function  without  poles  or  zeros  in  U2),:;sLJch 
that* |f*|  = |f*|.  (A  similar  statement  holds  for  Theorem  3).  However, 
the  proofs  of  these  converses  involve  some  technical  analytic  details,  and 
so  are  relegated  to  an  appendix. 

The  modification  in  the  above  converses  lies,  of  course,  in  the  fact 
thac  we  cannot  conclude  that  f itself  has  rational  spectral  factors;  thus 
there  are  some  rational  functions  which  can  be  stabilized  without  changing 
the  amplitude  response  but  to  which  the  classical  1-variable  factorization 
technique  cannot  be  applied.  A simple  example  of  this  is  the  function 

Zj  + Z2  - 1 

f(Zl’  Z2)  = Zj  + Z2  - ZjZ2 

Here,  |f*|  is  identically  1,  and  so  has  trivial  spectral  factors;  but  f it- 
self clearly  does  not. 

Although  the  converses  of  theorems  1 and  3 are  proved  in  the  appendix, 
there  is  another  result  related  to  the  converse  of  Theorem  3;  by 
strengthening  the  conditon  for  m=n=l  alone,  we  can  get  a stronger  converse 
for  polynomials.  Before  we  state  this  converse,  however,  we  first  give  a 
stability  criterion  (used  in  the  proof  of  the  converse)  which,  although 
previously  known,  [3],  has  not  appeared  in  the  engineering  literature. 

Although  not  as  sharp  (in  terms  of  dimension)  as  some  other  known  criteria  [4], 
it  has  two  advantages  which  make  it  useful  for  theoretical  purposes:  first, 
it  is  given  in  terms  of  a one-parameter  family  of  discs  without  the  lower- 
dimensional test  in  [5];  and  second,  unlike  most  other  stability  tests,  which 

conclude  the  nonvanishing  of  a polynomial  on  if2  from  its  nonvanishing  on 
__2  2 

some  subset  of  U which  contains  T , this  test  allows  the  polynomial  to 

2 

vanish  at  some  points  in  T , but  concludes  only  that  the  polynomial  does 

J 


» 
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2 

not  vanish  on  U . The  criterion  is: 

Theorem  4: 

Suppose  a polynomial  f has  no  zeros  in  the  set 

{(2,,  Z2)  t U2  | |Z, | = |Z2|); 

2 

then  f has  no  zeros  in  U . 

This  is  proved  in  a much  more  advanced  context  in  [3];  however,  it 
can  also  be  easily  proved  by  applying  one  of  the  criteria  in  [4]  to  the 
polydiscs 

= {(Zx.Zz)  e f2  | Ui|  1 r,  |Z2|  <_  r} 

for  o<r<l. 

For  the  hypotheses  imply  that  f has  no  zeros  on  the  distinguished  boundary 
? 

of  Ur  (for  o<r<l),  and  none  on  the  set 
{(Zj.Zz)  e (2|Z1=Z2>  n Vr 

? 

Thus  by  theorem  5 in  [4],  f has  no  zeros  in  Ur  for  any  r < 1,  and  so  f has 
2 

no  zeros  in  U . 

We  can  now  state  and  prove  the  partial  converse  to  theorem  3. 

Theorem  5: 

If  f is  a polynomial  with  the  property  that 

J log  f(e^0,e^0+f^)  jde  = log  |f(o,o)|.for  o<y<2v  , 

2 

then  f has  no  zeros  in  U . 

Thus  we  strengthen  the  condition  for  m=l  and  n=l  in  theorem  3 by 

specifying  that  the  constant  in  question  is  to  be  log  |f(o,o)|:  it  then 

follows  not  only  that  f has  rational  spectral  factors,  but  that  it  is 

2 

actually  zero-free  in  U . 


y 


By  theorem  4,  it  suffices  to  prove  that  f has  no  zeros  in  the  set 


{(Zr,Z2)  e U2|  IZJ  = |Z2|) 

But  this  set  is  the  union  of  the  open  discs 
UZi,Z2)|z2  = ej¥Zlt  | Zi | <1 1, for  o < » < 2*; 


we  therefore  wish  to  prove  that  f has  no  zeros  in  any  of  these  discs; 
or  equivalently  that  the  function  fy  of  one  variable  defined  by 

(Z)  « f(Z,ZeJ,p)  has  no  zeros  in  the  open  unit  disc.  Applying  Jensen's 
formula  [6,  p.299]  for  the  unit  disc  to  f^,  we  get 

h \Q  lo9  fT  (ej0)|de  = log  fy(o}|-J  log  )Z.| 

where  the  summation  is  over  all  the  zeros  (counted  with  multiplicity)  of 
fy  in  the  unit  disc.  Expressing  this  in  terms  of  f: 

2tt 

J0  lo9  j f(eJ0,eJ^0+y))  Jde  = log  f(o,o)|_  f log  |Zf| 
and  so  J]  log  |Z^|  = 0. 

Since  for  any  Z.  in  the  open  unit  disc  log  | Z . | < o,  the  conclusion  follows. 
(It  is  clear  from  the  proof  that  we  always  have 

(0  lo9  |f(eJ0,eJ^e+'l'))|de:  > log  |f(o,o)  j 

it  follows  from  this  that  in  fact  the  apparently  weaker  condition 
2ir  ^2f 

4TT2  lo  lo  109  f(eJ01  »eJ02)|d0id02  = log  |f(0.0)| 

2 

is  sufficient  to  guarantee  that  f is  zero-free  in  U . See  [2,p.73]). 

Stable  IIR  Filters  and  Minimum-phase  FIR  Filters: 

The  very  close  relationship  of  spectral  factorization  to  the  nonvanishing 
2 

of  polynomials  in  U , and  thereby  to  stable  IIR  filters  (via  the  denominator 
polynomial)  and  minimum-phase  filters  (via  the  numerator  polynomial)! is  al- 
ready clear  from  the  previous  sections.  The  force  of  theorem  2 is  that 
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purely  from  the  point  of  view  of  amplitude  response, transfer  functions  having 

rational  spectral  factors  are  equivalent  to  those  without  poles  or  zeros  in 
2 

U . Thus  the  restrictions  on  amplitude  response  in  theorems  1 and  3 apply 
to  the  denominator  polynomial  of  any  stable  1 1 R filter;  the  contribution  of 
the  denominator  polynomial  to  the  overall  amplitude  response  of  the  filter 
(in  the  case  of  an  all-pole  filter,  the  entire  amplitude  response)  must 
satisfy  the  restrictions  imposed  by  theorems  1 and  3.  We  have,  therefore, 
identified  the  properties  of  the  amplitude  response  which  make  it  impossible 
to  stabilize  a filter;  if  the  original  amplitude  response  has  large  ,'overal  1 
variation  in  the  "wrong"  directions,  attempting  to  find  a stable  filter  which 
closely  matches  this  response  is  futile.  Close  matching  of  the  amplitude 
forces  instability.  This  has  already  been  shown  by  example  by  Bose  [9]  and 
Woods  [10]:  we  now  see  that  it. is  the  variations  in  the  amplitude  response 
in  the  "wrong"  directions  in  their  examples  which  accounts  for  their  be- 
haviour. 

It  is  also  of  interest  to  note  that,  in  the  Shanks  procedure  of  minimizing 
||  |fg-l|2  dexdQz 

over  all  polynomials  f of  given  degree  (where  g is  the  original  polynomial), 

if  the  allowable  f’s  were  restricted  to  those  which  have  polynomial  spectral 

2 

factorizations,  the  procedure  would  yield  a polynomial  devoid  of  zeros  in  U . 

It  does  not  appear  that  this  observation  can  be  used  as  the  basis  for  a work- 
able stabilization  method,  however,  since  the  condition  that  f have  polynomial 
spectral  factors  is  intractabl y non! inear  in  the  coefficients  of  f;  and  further, 
in  many  cases  this  procedure  would  yield  an  f which  was  only  marginally  stable. 
For  the  same  reasons,  restricting  oneself  throughout  the  design  procedure  to 
polynomials  which  satisfy  the  condition  in  theorem  5 does  not  appear  to  be  a 
feasible  method  of  ensuring  stability. 
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Examples  and  Comments: 

An  example  of  the  behaviour  of  those  polynomials  not  possessing  poly- 
nomial spectral  factors  has  already  appeared  in  the  literature,  although 
in  a different  context;  we  repeat  this  example  here. 

A(Z!,Z2)  = 1 75ZX  + .9Zi  + 1.5Z2  - 1.2ZxZ2  + 1 .3zfz2  + 1.2Z2  + .91^1  + .5ZxZ2 

This  polynomial  was  studied  in  [7];  the  associated  Shanks  polynomial  was 
found  to  be  stable  but  to  have  a substantially  different  amplitude  response  • 
from  that  of  A (for  more  details,  see  [7]).  The  fact  that  A does  -not  have 
polynomial  spectral  factors  was  established  by  checking  the  condition  in 
theorem  3 for  m=n=l  and  f = o,  V * ir,  with  the  following  results:  (correct  to 
nine  decimals) 

^ L"109  l^0’  ej9)|de  = .696570700 

2tt 

jQ  log  |A(eJ0,  ej(6+^)  de  = 1.134686936 

As  an  example  of  a polynomial  with  rational  spectral  factors,  we  have 

B(ZX  ,Z2)  = 1 + 2.25ZX  + 2.25Zz  + .5Z*  + .5Z2  - 6.51^2  - Z^Z2  - ZXZ2 

2 2 
- 4ZiZ2 

This  factors  into  (1  + .25Zx  + .25Z2  + .5ZXZ2)(1  + 21x  + 2Z2  - 8ZjZ2) 

2 2 

the  first  factor  factor  having  no  zeros  in  U , the  second  none  in  V ; re- 

2 

versing  the  second  factor  gives  a polynomial  without  zeros  in  U : 

B(Zi,Z2)  = (1  + .25ZX  + .25Z2  + .5  ZxZ2) (-8  + 2Z2  + 2Zx  + ZXZ2) 

2 2 2 2 2 2 
= -8^-2ZiZ2  + .5  Zx  + .5Z2  + 1.25ZxZ2  + 1.25ZXZ2  + .5ZXZ2, 

*\» 

and  B has  the  same  amplitude  response  as  B. 
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In  order  to  gain  some  idea  of  the  stringency  of  the  conditions  in 
theorem  3,  let  us  consider  the  case  of  an  ideal  band-pass  filter.  By 
an  ideal  band-pass  filter  we  will  mean  a filter  whose  amplitude  response 
is  equal  to  1 on  some  subset.  A,  of  the  square  o<91<2ir,  o<e2<2ir,  and 
equal  to  K « 1 on  the  complement  of  A (of  course  this  specification  con- 
tinues over  the  whole  plane  by  periodicity).  This  of  course  is  not  the 
amplitude  response  of  any  rational  function,  but  in  practice  for  certain 
shapes  of  the  set  A,  one  may  wish  to  approximate  such  a response  by  a 
rational  function.  One  easily  sees  that  up  to  a scale  factor,  the  averages 
in  theorem  3 are  in  this  case  merely  the  fraction: 

length  of  the  line  L.  lying  in  the  complement  of  A 
Total  length  of  the  line  L. 

It  is  easily  seen  from  this  that  there  are  very  few  passband  shapes  of 
practical  interest  which  satisfy  even  the  first  of  these  conditions  (where 
n=l  and  m=l);  in  other  words,  very  few  which  can  be  accurately  approxi- 
mated by  transfer  functions  having  rational  spectral  factors.  (This  is 
not  to  imply  that  one  would  in  practice  be  restricted  to  such  filters;  the 
above  discussion  is  meant  solely  as  an  indication  of  the  severity  of  the 
restrictions  on  the  amplitude  of  such  filters). 

Finally,  we  remark  that  there  does  not  seem  to  be  any  difficulty  in 


I 


extending  the  results  in  this  paper  to  higher  dimensions,  and  to  multi- 
dimensional systems  other  than  digital  filters. 


72 


APPENDIX 


The  converses  to  Theorems  1 and  3. 

These  converses  involve  some  technical  ideas  and  results  from  [2];  the 
most  important  ideas  are  those  of  inner  function  [2 , p.  1 05] , outer  function 
[ 2,p.  72],  Poisson  integral  [2,  p.  1 7 ] and  the  classes  N(U2)  [2,  p.44]  and 
MU2)  [2,  p.44]. 

We  will  also  use  the  following  notation  from  [2]  (Here  f is  an  analytic 
function  on  U ): 

i.  f*(eJ01,  e^02)  4 lim_  f(reJ01,  re^02  ) 

r-*-l~ 

will  denote  the  radial  limit  of  f 

(this  is  clearly  consistent  with  our  previous  use  of  f*); 

2 

ii.  For  w = (w^  w2)  e T , fw(Z)  will  denote  the  one-varible 
function  defined  by 

fw(Z)  4 f(Zwlf  Zw2) ; 

2 

iii.  if  (])  is  a function  defined  on  T which  is  absolutely  integrable 

there, 

: 1 

4*(ni»n)4  2 2 exp(.jm6l  -jne2)  (J>  (e i , e 2 ) de2de1 

. o o 

will  denote  the  Fourier  coefficients  of  ({). 

2 

iv.  For  any  function  <p  on  T , 

• * 2 IT  C 2 TT 

7T7  <t>(0i»02)  d02d0i  w-*^  be  denoted  by 

4ir  o JO 


$ dm  or  <()  (w)  d m (w) 

JT2  JT2 


We  will  first  prove  the  converse  to  Theorem  1,  and  from  this  derive 

i 

the  converse  to  Theorem  3.  First  of  all,  however,  we  need  the  following 


lemma  (which  is  given  as  a problem  in  [2]). 


Lemma  A1 : 


If  $ is  a real -valued  function  defined  on  T2 
such  that 

t e L (T2)  (i .e. , ( | (p  |dm  < «) 

Jj2 

and 


<j)(m,n)  = 0 for  mn  < 0, 

then  there  is  an  outer  function  f on  U2  such  that 
P[<1>]  = log  | f | 

(where  P[  ] denotes  "Poisson  integral  of"). 


Proof 


Let  a 


mn 


and  let 


$(m,n) 

1/2  <j)(m,n) 

OO  00 


(m,n)  ^ (o,o) 
(m,n)  = (o ,o ) 


9(Z,.Z2)  -i  E amn  Z5 
m=o  n=o 
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This  series  clearly  converges  uniformly  on  compact  subsets  of  U2, 
and  so  defines  an  analytic  function  there. 


/ 
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If  we  let  f = e9 


then  f is  analytic  in  U , and 


log  | f | s I z amn  rf  r2  exp  (jme:  + jne2) 


m=o  n=o 


+ 2 1 arnn  ri  r 2 exP  (-jmfli  -jne2) 

m=o  n=o 


* E E $ (m,n)rjmMnUxp(jme1  + jn02  ) 
m=-<»  n=-« 

- P[<P]  [2„p.l7] 

Next  we  prove  that  f is  outer;  we  have  (for  0<r<l) 


| 2log+|f(rw)|dm  (w) 


< 

< oo 


t2  1 1 °g  I f (r»w)l  |dm ( w ) 
^2  I pC(J)3  ( rw)  | dm ( w ) 

|^2|$(w) )dm(w) 


[2 ,Thm. 


and  so  f e N(ir). 

Now  f*  exists  almost  everywhere  on  T2  [2,Thm.3.3.5]  and  log|f*|  = 
everywhere  on  T2  [2,Thm.2.2.1];  thus  log  |f|  = P[log  |f*|]  and  so 
f e N*(U^)  [2,Thm.3.3.5],  and  log  |f(0)|  * J 2 log  | f*(w) |dm(w) . 

Thus  f is  outer. 

Q.E.D. 


2.1.3(c)] 


<|)  almost 


We  can  now  prove  the  converse  to  Theorem  1 : 


» 
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Theorem  A2: 


Let  ffZj.Zg)  be  a rational  function  (£0),  and  let 
<f>  = log  |f*| 

A 

If  0(m,n)  * 0 for  mn<0,  then  there  is  a rational  function  g without  poles 
or  zeros  in  such  that  |g*|  * |f*|. 

Proof; 

By  Lemma  Al,  there  is  an  outer  function  g such  that 
log  ] g | = P[log|f|]. 

This  implies 

p 

log  |g*|  = log  |f*|  almost  everywhere  on  T . 

p 

Therefore,  for  almost  all  wcr 

log  |g*  (Z)|  * log  |f*  (Z)|  for  almost  all  ZeT  [2, Lemma  3.3.2], 

p 

and  gw  is  outer  for  almost  all  w eT  [2, Lemma  4.4.4]. 

For  any  such  w,  let  Zj,  ...,  Zn  denote  the  poles,  and  Zn+^ , . ..  Zm  the 
zeros,  of  fw(Z)  in  U,  and  let 

' n 7-7  m 7 7-1 

fw(z)  = n iik  n zk  1 1 f (z) 

" i i i ■ y i V» 


k=l 


2.  z-i 


k=n+l  Z-lL 


Then  f has  no  poles  or  zeros  in  U and  is  rational;  hence,  f is  outer.  Since 

n W 

gw  is  outer,  we  have  fw/gw  is  outer.  Also  |fw*|  - |fw*|,  and  so  |fw*|  * ^*1 
for  almost  all  ZeT.  Thus  fw/9w  is  inner.  But  a function  which  is  both  outer 
and  inner  is  a constant  of  modulus  1,  and  so 


’w 


w 


for  some  real  f. 
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f 

' 

! 


o 

Thus  gw  is  rational  for  almost  all  weT  , and  so  gw  is  rational  for  all 

w e E,  where  EcT  is  a compact  set  of  positive  measure  (by  the  inner 

regularity  of  the  measure).  It  follows  by  [2,Thm.  5.2.2]  that  g is 

rational  (since  the  vanishing  of  a polynomial  P on  a set  of  positive  measure 
2 

in  T would  imply 

log  |P*|  i L](T2) 

and  so  P = 0.) 

2 

Thus  g is  a rational  function  without  poles  or  zeros  in  U , and 

| g*|  * | f* | almost  everywhere  in  T2 
and  so,  since  g and  f are  both  rational, 

|g*|  ■ |f*|  on  T2. 

Q.E.D. 

We  next  prove  the  converse  to  Theorem  3: 

Theorem  A3: 

Let  f(Z-|  .Zg)  be  a rational  function  0)  and  let 


If 


1 f2* 

2-  L 


(p  = log  |f*| 


(J)  (me,  ne  + f)de  is  a constant  independent  of  v for  each 


pair  (m,n)  with  m >o  and  n>o  then  there  is  a rational  function  g without 
poles  or  zeros  in  U2  such  that  |g*|  = |f*|. 


i 


* 
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Proof: 


Let  m > o,  n > o,  and  let  J.  f o be  an  integer. 


Then 


rZTT 

; jimf  $(me,ne +'?)de  dv  = 0 
J o e Jo 

1 21*  (2i 1 . 

Jo  Jo  eJim41  $ (me.ne +’p)d8dH'  * 0 


Making  the  change  of  variables  defined  by 

9 = 1 0 1 
m 

y = 92-  n_  0J , 


m 


we  get 


1 f21T  fffl01+2lT 

m Jo  JlL*  exP  (jin>02-  itneO  $ (0i ,Q2)de2^i 

m 1 

and  since  the  integrand  is  periodic  in  0!  and  e2 

(2*  [2  it 

J exP  (jJlm02.  - jjJn© i ) $ (0i,02)de2d0i  - 0 
o o 

and  so  (jj  (-in,  im)  = 0 for  all  Ji  f 0,  m > o and  n > 0, 
that  is, 

A 

(Jl(m,n)  = 0 for  all  m,n  with  mn  < o. 


* 0 


The  result  now  follows  from  Theorem  A2. 


Q.E.D. 


Finally,  we  note  that  if  f in  Theorem  A3  is  a polynomial,  then  the  converse 


In  Theorem  2 Implies  that  f has  polynomial  spectral  factors.  Thus  we  have  the 
full  converse  of  Theorem  3 for  polynomials. 
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SYMMETRIC  HALF-PLANE  FILTERS 


John  Murray 


ABSTRACT 

A class  of  two-dimensional  recursive  digital  filters  called  symnetric 
half-plane  filters  is  discussed;  some  properties  of  these  filters  are  derived 
and  it  is  shown  that  in  certain  situations  these  properties  may  give  the 
symnetric  half-plane  filters  both  theoretical  and  practical  advantages  over 
previously  proposed  filters.  In  particular,  they  are  ideally  suited  to  highly 
parallel  processing. 

INTRODUCTION 

In  the  literature  on  2-dimensional  recursive  digital  filters,  two  main 
types  of  filter  have  been  studied;  these  are  the  quarter-plane  filter  (e.g.[l], 
[2])  and  the  asymmetric  half-plane  filter  [3].  Basically,  the  two  correspond 
to  two  different  concepts  of  causality.  The  general  stability  conditions 
for  a wide  class  of  filters  (including  symmetric  half-plane)  were  discussed  in 
[4];  unfortunately,  however,  those  filters  are  not  recursively  implementable  in 
general.  Here  we  will  consider  a class  of  filters  which  are  recursively  im- 
plementable, and  satisfy  the  same  stability  conditions  as  those  in  [4]. 

SYMMETRIC  HALF-PLANE  FILTERS 


By  a symnetric  half-plane  filter  we  will  mean  a (causal,  recursive)  2-di 
wenslonal  digital  filter,  the  denominator  of  whose  transfer  function  is  of 


t 
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This  differs  from  the  filters  in  [4]  in  that  m goes  from  1,  rather  than  0, 
to  M;  i.e.,  this  filter  omits  all  of  the  row  m=0  except  for  the  constant 
term;  the  asymmetric  half-plane  filters  omit  half  of  this  row.  The  filter  (1) 
is  recursively  realizable,  since  the  computation  of  the  output  at  any  point 
depends  only  on  the  outputs  in  previously  computed  rows;  looked  at  from  another 
point  of  view,  each  row  of  output  depends  only  on  previous  rows  of  output.  This 
has  two  effects;  firstly,  it  focuses  attention  on  the  row  as  the  basic  element 
in  the  filter;  secondly,  it  implies  that  all  the  outputs  in  a given  row  may  be 

computed  in  parallel,  since  each  output  in  a row  depends  only  on  outputs  in 

previous  rows,  and  not  on  any  of  the  outputs  of  the  same  row.  This  is  the  main 

practical  advantage  of  this  class  of  filters  - it  would  be  of  significance,  how- 

ever, only  in  real-time  hardware  applications  of  2-dimensional  filtering,  and 
these  seem  to  be  few. 

SOME  PROPERTIES 

Using  the  methods  in  [4],  one  can  easily  derive  the  following: 

The  filter  (1)  (i.e.,  the  all-pole  filter  whose  denominator  is  A(Zj,Z2))  is 
stable  if  A(Z],Z2)  / 0 for  all  (Z]t  Z2)  such  that  |Z] | =1  and  |Z2|  <1. 

We  note  that  this  set  is  the  same  as  that  for  the  symmetric  half-plane  filter 
in  [4];  it  is  smaller  than  that  for  the  asymnetric  half-plane  filters  [3].  It 
is  the  smallest  "instability  set"  (known  to  the  author)  of  any  recursively  im- 
plementable  class  of  filters. 

However,  there  is  a price  to  be  paid;  the  amplitude  response  of  the  filter 
is  restricted  as  follows: 

If  A(ZpZ2)  is  of  the  form  (1),  and  if  ^ ^ is  the  transfer  function  of  a 

stable  filter,  then 

log  |A(e^01 ,e^02) | de2  = 0 (2) 

■<0 


independently  of  0i. 
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Thus,  the  average  gain  along  any  line  of  length  2ir  parallel  to  the  e2-axis  is 

constant;  or  in  other  words  the  filter  cannot  have  variations  in  the  ej-direction 

in  overall  (average)  gain.  Equivalently,  if  the  cepstrum  of  |A(eJ'®1,  eJ02)|  is 

given  by  l l amn  Zl  z*  , 
m=-«  n= *-• 

<v 

then  a„„  = 0,  for  all  n. 
on 

This  follows  immediately  from  (2)  and  the  definition  of  the  cepstrum;  (2) 
will  be  proved  in  a forthcoming  paper. 

This  implies  that  in  order  to  realize  an  arbitrary  magnitude  function,  the 
filter  must  either  have  a (nonminimum-phase)  numerator  polynomial,  or  the  filter 
must  be  cascaded  with  a 1-dimensional  filter  in  Z-j . It  is  very  easy  to  calculate 
the  ideal  amplitude  response  of  this  filter. 

DESIGN  AND  IMPLEMENTATION  CONSIDERATIONS 

It  is  conceptually  convenient  (and  in  a large  number  of  cases,  computation- 
ally efficient)  to  implement  the  Convolution  in  the  Z^-direction  by  means  of  the 
Fourier  Transform.  (It  is  assumed  from  the  beginning  that  the  dimension  of 
the  array  to  be  filtered  is  a known  fixed  constant  in  the  Z^ -direction,  i.e.,  each 
row  is  of  the  same  fixed  width).  From  this  point  of  view,  and  regarding  each  row 
as  a single  entity  described  by  its  1-dimensional  Z-transform,  the  coefficients 
dmn  ‘n  are  irrelevant;  what  matters  are  the  M functions 

V*j6l>  “ 1 a-(eJ9,,N- 

Further,  the  stability  requirement  for  the  filter  is  equivalent  to  the  requirement 
that  for  each  fixed  elf  the  filter  defined  by 

] 

E Sm(ej9>)Z?  ♦ 1 

m=l 

Finally,  the  functions  a (eJQl)  do  not  have  to  be 


is  1-dimensionally  stable. 
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analytic  or  meromorphic  functions;  this  is  seen  by  letting  N-h».  In  other  words, 
the  roots  of  am(eJ 0 1 ) can  vary  quite  arbitrarily  with  ej.  Thus,  we  can  design 
the  one-variable  filter  A (Z2)  by  any  one  of  the  usual  one-variable  design 
methods  we  choose  (yielding  a stable  filter)  for  each  Qj : the  result  will  be  a 
stable  two-variable  filter:  further,  if  our  one-variable  method  gives  poles  and 
zeros  explicitly,  we  have  the  same  for  our  two  variable  filter,  which  can  there- 
fore be  expressed  as  a cascade  of  filters  of  degree  1 in  Z2.  Finally,  if  one 
desires  a filter  of  finite  degree  in  Zj,  one  can  solve  the  following  approximation 

problem  (for  each  m,  1 < m < M);  minimize  (over  bl 

— — mn 

n 

II  E bmnejn91  " t>n/ej9lHI  subject  to 
n=-N 

| y b jne x | < 1 for  all  elt  where 
n--N  mn 

bm(e^01)  denotes  the  m-th  pole  of  A0 ^ ( Z2 ) as  a function  of  elf  and  ||  ||  denotes 
some  error  norm. 

Hopefully  this  will  become  clearer  on.  consideration  of  the  following  example. 
EXAMPLE 

We  wish  to  design  a filter  with  second-order  Butterworth  response  in  Z2  to 
approximate  the  fan  filter  whose  passband  is  the  set  |e2|  < 1 8 1 1 . For  fixed  elt 
therefore,  the  filter  is  a 1-dimensional  filter  whose  passband  is  the  set  |e|  < lej 
Using  the  bilinear  transform  technique,  we  find  the  second-order  continuous 
Butterworth  filter 

1 

1 + /2  s/ta>c  + S2/u>2 


transforms  into 
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"c*  (1  ♦ Zj)1 

*c  **e  * 1 ♦ (*»e*  - 0 Zj  ♦ (»c*  ♦ /J  »e  ♦ 1)  zj 


(3) 


while  its  (stable)  poles  transform  into 

11  + -c'TJ11'!/'!  - •c'tH’ 

and 

tor  m < > 0. 

c — 

A 

In  accordance  with  the  usual  frequency  warping,  we  take  u>c  3 |tan  1/2|;  however, 
we  note  that  this  causes  stability  problems  at  0i=O  and  ex *ir ; we  therefore  take 
a small  perturbation  ofu>c,  e.g.. 


_C 

C 


(4) 


The  filter  can  now  be  directly  implemented  by  multiplying  the  Fourier  transform 
of  the  previous  output  rows  by  the  appropriate  functions  according  to  (3)  and  (4) 
and  performing  the  recursion  from  row  to  row  directly. 

Alternatively,  if  a finite-order  (in  Zt)  filter  is  desired,  we  must  solve 
the  problem:  Minimize  (over  bln) 


*V;i  - ? 

* */2 1 J-,  " 


subject 


to  ' " I * 1 for  a11 

n*-N 


and  similarly  for  the  other  root.  The  resulting  filter  may  then  be  implemented  in 
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cascade.  It  should  be  pointed  out  that  while  the  above  optimization  problem  is 
not  simple,  it  is  one-dimensional. 

Finally,  the  above  example  was  chosen  for  simplicity  and  convenience  rather 
than  realism.  Clearly  a Butterworth  filter  is  not  optimal  for  this  problem, 
especially  when  it  yields  a design  which  is  not  all-pole;  if  we  have  to  store 
input  rows  (as  well  as  output  rows)  we  may  as  well  use  them,  and  design  an 
elliptic  type  filter;  alternatively,  we  might  use  a filter  which  is  all-pole 
in  its  discrete  form. 
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Introduction  -> 

The  concepts  delineated  in  this  paper  arose  in  part  from  an 
introductory  study  of  Riemann  surfaces.  Associated  with  an  analytic 
function  is  a Riemann  surface.  It  has  the  property  that  the  image 
of  simply  connected  regions  in  the  complex  plane  are  simply  connected 
on  the  Riemann  surface. 

The  point  made  here  is  that  the  Nyquist  criterion  is  trivial 
for  simply  connected  regions.  Moreover,  if  one  can  work  on  the 
Riemann  surface,  this  trivality  carries  over  to  the  general  case. 

To  illustrate  the  point,  let  Figure  1-a  be  the  image  of  the  right 
half  plane  under  an  analytic  map.  The  region  is  not  simply  connected. 
Figure  1-b  shows  the  "same  region"  as  it  might  appear  on  an 
appropriate  Riemann  surface.  Here  the  region  is  simply  connected. 

t 


I 


Figure  1 


Under  the  hypothesis  that  f is  bounded  at  infinity,  the  boundary 


of  the  region  in  Figure  1-b  is  the  image  of  the  imaginary  axis  as 
indicated  by  the  darkened  line  in  Figure  1-a. 
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Now  remove  "-1"  (this  may  be  a set  of  points)  from  the  Riemann 
surface.  The  essential  argument  we  need  is  that  the  Nyquist  contour 
in  the  complex  plane  is  homotopic  to  zero  if  and  only  if  "-1"  is  in 
the  interior  of  its  image  on  the  Rimann  surface. 

Although  motivated  by  the  intrinsic  properties  of  Riemann  surface, 
this  paper  drops  any  further  discussion  of  the  concept  so  as  to 

simplify  the  exposition.  Instead,  the  paper  exploits  the  fact  that 
the  Nyquist  contour  is  a simple  closed  curve  in  the  complex  plane. 
Mathematically  we  draw  only  on  the  intuitive  concept  of  homotopic 
triviality  as  found  in  algebraic  topology. 

In  the  sequel,  we  prove  the  classical  stability  results  via 
homotopy  theory.  In  particular,  we  utilize  covering  space  theory. 

We  believe  our  analysis  is  clearer  and  more  intuitive  than  has 
hitherto  appeared.  Moreover,  we  believe  that  this  research  indicates 
that  the  nub  of  the  Nyquist  criteria  is  in  fact  homotopy  theory. 

In  a future  paper,  we  will  generalize  these  results  to  functions  of 
several  complex  variables  and  their  application  to  the  stability  of 
of  multi-dimensional  digital  filters. 

/ 

Mathematical  Preliminaries 

O 

Firstly,  let  £,  £+,  and  £+  be  the  complex  plane,  the  closed 
right  half  plane,  and  the  open  right  half  plane  respectively.  Let 

O 

~ l - £+.  Basic  to  homotopy  theory  is  the  concept  of  a path. 

A path  or  a curve  in  the  complex  plane  is  a continuous  function 
Y*-[0 »!]-*£.  For  this  paper  all  paths  are  understood  to  be 
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rectifiable— i.e.  y is  also  a function  of  bounded  variation.  (2) 
y is  a closed  path  if  y(0)  3 yO)*  y is  a simple  closed  path  if  y 
is  a closed  path  and  has  no  self  intersections.  The  image  of 
I 3 [0,1]  under  y is  called  the  trace  of  y and  is  denoted  by  (y>- 

Two  closed  curves  yq  and  y-j  are  homotopic  (y0  - y-|)  in  Z if 
there  exists  a continuous  function  r:  I x I jt  such  that: 

(a)  r(s,0)  - y0(s)  0 < s < 1 

(b)  r(s,l)  3 yi(s)  0 _<  s <_  1 

(c)  r(0,t)  =*  r(l,t)  0<t<l 

Intuitively,  yq  is  homotopic  to  Yj  if  one  can  continuously 
deform  yq  into  Y] • Moreover,  it  is  easily  shown  that  the  homotopy 
relation  is  an  equivalence  relation.  (4)  (5) 

Another  important  property  of  a closed  curve  is  its  index.  The 
index  of  a closed  curve,  y,  with  respect  to  a point  "a"  not  in  {yl 
Is: 

n(Y'*a)  = / (z-a)  ^dz  (2) 

Y 

Observe  that 

/(z-aj^dz  = Jd(ln(z-a))  = Jd(ln|  (z-a)  | ) + i/d(arg(z-a) ) 

Y Y Y Y 

= i/d(arg(z-a)) 

Y 

This  integral  therefore  measures  i times  the  net  increase  in  angle  that 
the  ray  r of  Figure 3 accumulates  as  its  tip  traverses  the  curve  y. 
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Figure  2 


Following  the  comments  of  J.  Barman  and  J.  Katzenelson  (1), 
for  the  integral  to  be  well  defined  it  is  necessary  to  specify  the 
appropriate  branch  of  arg(z-a)  at  each  point  of  the  integration. 

We  will  assume  the  choice  of  branch  as  outlined  in  (1). 

Finally,  we  point  out  that  this  definition  of  index  (encirclement) 
is  a special  case  (i.e.  in  the  plane)  of  the  general  topological 
concept  of  Brouwer  degree.  (4)  (5)  (6) 

At  any  rate  n(\;a)  = 0 if  and  only  if  y is  homotopic  to  a point 
in  l - {a},  (cf.  prop  5.4,  ref.  2)  Simply  then,  a closed  curve  y 
does  not  encircle  the  point  "-1"  if  and  only  if  y is  homotopic  to  a 
point  in  t - {-1>.  We  will  henceforth  refer  to  such  a y as  being 
homotopically  trivial. 

Conversely,  y encircles  "-1"  if  and  only  if  y cannot  be  continuously 
deformed  to  a point  in  l - {-1}.  Clearly  these  ideas  indicate  that 
the  Nyquist  encirclement  condition  is  fundamentally  a homotopy 
concept. 

To  further  illucidate  the  point,  let  g(s)  be  a rational  transfer 
function  depicting  the  open  loop  gain  of  a scalar  single  loop  feedback 

A O A 

system.  Suppose  all  poles  of  g(s)  are  in  and  g(»)  < m < «.  Via 
the  Nyquist  Criteria,  the  closed  loop  system  is  stable  if  and  only 
if  h(s)  = g(s)/(l+g(s))  is  stable;  if  and  only  if  the  image  of  the 
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A A 

imaginary  axis,  under  g(s),  (the  Nyquist  plot  of  g(s) ) does  not  pass 
nor  encircle  "-1". 

Specifically,  the  encirclement  of  "-1"  by  the  Nyquist  plot 

o * 

implies  there  exists  at  least  one  Sg  in  0+  such  that  g(Sg)  = -1. 

Thus  the  Nyquist  contour  is  homotopically  trivial  in  0+  - {<T^(-1)} 
if  and  only  if  the  Nyquist  plot  is  homotopically  trivial  in 

•gte+M-i). 

Motivation  for  this  approach  also  arose  from  a close  scrutiny 
of  the  classical  proof  of  the  Nyquist  criteria  which  depends  on  the 
argument  principle.  The  argument  principle  supplies  unnecessary 
although  specific  information  in  that  it  counts  the  number  of  times 
"-1"  is  encircled.  This  may  account  for  the  apparent  difficulty  in 
generalizing  the  Nyquist  criteria.  Nevertheless,  the  affinity 
between  homotopy  and  encirclement  ideas  led  the  authors  to  a minor 
study  of  algebraic  topology. 

In  our  setting,  algebraic  topology  establishes  a topologically 
invariant  relationship  between  a metric  space,  X,  and  an  algebraic 
group  called  the  fundamental  group  of  X,  denoted  by  ir(X).  The 
relationship  is  topologically  invariant  in  that  homeomorphic  spaces 
have  isomorphic  fundamental  groups. 

Specifically,  the  fundamental  group  is  a set  of  equivalence 
classes  of  closed  curves.  Each  equivalence  class  consists  of  a set 
of  curves  homotopically  equivalent.  The  group  operation  is 
"concatenation"  of  curves. 


t 
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For  example,  the  fundamental  group  of  £ consists  of  one 
element,  i^,  the  identity,  since  all  closed  curves  are  homotopic  to 
zero.  If  X = t - (-1).  then  ir(X)  has  a countable  number  of  elements: 
ix  (the  identity)  equal  to  the  equivalence  class  of  all  closed  curves 
not  encircling  "-1"  and  the  remaining  elements,  un  (n  * 1,  2,  3 ...) 
consisting  of  the  equivalence  class  of  all  closed  curves  encircling 
"-1",  n times.  Moreover,  u-  concatenated  with  is  equal  to  the 
element  . 

Now  let  X and  Y be  metric  spaces.  Let  f:  X -*■  Y be  locally 
homeomorphic.  In  particular,  assume  that  for  each  point  y in  Y there 
exists  an  open  neighborhood  G of  y such  that  each  connected  component 
of  f“^(G)  is  homeomorphic  to  G under  the  map  f.  Under  this  condition 
X is  said  to  be  a covering  space  of  Y.  (2)  (4)  Also  let  tt(X)  and 
it ( Y ) be  the  fundamental  groups  associated  with  X and  Y respectively. 
With  these  assumptions,  f effects  a group  isomorphism  (i.e.  a one  to 
one  onto  mapping  preservinggroup  operations)  <pf  between  tt(X)  and  a 
subgroup  of  tt(Y)  as  in  the  following  diagram.  (4)  (5) 


Figure  3 


96 


F is  the  functor  which  establishes  the  relationship  between  a 

topological  space  and  its  fundamental  group. 

Before  judiciously  tailoring  the  complex  plane  so  as  to  apply 

the  above  result,  we  distinguish  between  a critical  point  and  a 

critical  value.  A point  zQ  in  Z is  a critical  point  of  a differentiable 

function  f if  f'(z0)  =0.  A critical  value  of  f is  any  point  w = 

f(zQ)  whenever  zQ  is  a critical  point. 

Now  suppose  f:Z-*-Z  is  a rational  function  whose  set  of  poles 

is  P = {P-|,  ...»  Pn>.  Let  Q = {q-j q^}  be  the  set  of  all  points 

in  1 such  that  f(q^)  is  a critical  value  of  f.  Note  that  there  may 

be  q^'s  which  are  not  critical  points.  To  see  this  consider  g(z)*  = 

2 

z (z-a).  gtO)  = 0 implies  "0"  is  a critical  value  of  g,  but  g(a)  = 

0 with  g'(a)  t 0.  Finally,  define  T = {t^  = f -1  (-1 ) , i=l,...,n>. 

Note  also  that  since  f is  a rational  function,  P,  Q and  T are  finite 
sets.  Define  X = Z - {P  Q T}  and  define  Y = f(X). 


Lemma  1:  Under  the  above  hypothesis,  X is  a covering  space  of  Y. 
Proof:  For  X to  be  a covering  space  of  Y,  each  y in  Y must  have  an 
open  neighborhood  Gy  such  that  each  component  of  f"^(Gy)  is  homeo- 
morphic  to  Gy.  Using  the  inverse  function  theorem  (7  ) we  construct 
such  a neighborhood. 

Let  { x i xk>  = f ^(y)  where  again  the  finiteness  of  this 

set  is  a consequence  of  the  rationality  of  f.  Let  W-j , ...,  Wk  be 
disjoint  open  neighborhoods  of  x-j , ...,  xk  respectively.  Since 
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f is  analytic  on  X and  since  f‘(x)  f 0 for  all  x in  X,  the  inverse 

function  theorem  guarantees  that  there  exist  open  neighborhoods 

■ 

U-jCWj  (i  * 1»  .-..k)  such  that  is  homomorphic  to  * f (U - ) » 
where  it  follows  that  V1-  is  an  open  neighborhood  of  y. 

Thus  VR)  = Un.  Define  V = V]  fY...f*Vk. 

Clearly  V is  an  open  neighborhood  of  y and  f^OOcUjU  ...  U Uk- 
Since  each  is  homeomorphicto  f_1(V)C.Uj  is  homeomorphic 

to  V. 

Therefore  each  y in  Y has  an  open  neighborhood  Gy  such  that  * 
f (Gy)  has  each  of  its  components  homeomorphic  to  Gy.  It  follows 
that  X is  a covering  space  of  Y. 

Corollary:  The  fundamental  group  ir(X)  of  X is  isomorphic  to  a sub- 
broup  N of  ir(Y) . 

This  corollary  says  that  a closed  curve  in  X is  homotopically 
trivial  if  and  only  if  its  image  under  f is  homotopically  trivial. 

The  Scalar  Nyquist  Criterion 

In  this  section  we  apply  the  above  corollary  to  an  "ugly" 

Nyquist  contour.  After  proving  the  Nyquist  Theorem  using  this  "ugly" 
contour  we  relate  it  to  the  usual  Nyquist  contour.  This  will 
establish  the  classical  result. 

A 

Let  g(s)  be  a rational  function  which  represents  the  open  loop 

gain  of  a scalar,  single-loop  unity  feedback  system.  We  assume 

A 

g(s)  t 0.  Thus  the  closed  loop  system  has  a transfer  function 


— - ■*  .- ... - - 
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h(s)  » g(s)(l+g(s))-1. 

We  will  say  that  the  closed  loop  system  h(s)  is  stable  if 

a 

and  only  if  h(s)  has  all  its  poles  in  and  h(-)  < - . 

A 

Let  P = {p-j,  ....  pn>  be  the  set  of  poles  of  g(s)  and  let 
Q = (q-j,  •••»<!„,}  be  the  set  of  points  q^  such  that  g(q^ ) is  a critical 
value  of  g.  Define  T = {t^ | t^  = f-1(-l),  i = 1,  ...»  1).  Finally, 
let  X = £ - {POQuT}  and  let  Y * f(X).  Lenina  1 implies  X is  a 
covering  space  of  Y under  the  mapping  f. 

A 

Assume  for  the  present  that  g(iw)j*  "-1"  for  -•£«<-  . The 


first  task  is  to  construct  the  "ugly"  Nyquist  contour  as  well  as 
the  classical  contour.  Define  the  ulgy  contour  to  be  yR  where 
yr*.  I -*•  Xc£  is  a path  whose  trace  is  illustrated  in  Figure  4-2. 

Note  that  R is  chosen  strictly  greater  than  max  (|p^|,  [q^.  |,  |tk() 
for  1 <.  i ± n,  1 <.  j <.  m,  and  1 £ k £ 1.  The  indentations,  along 
the  imaginary  axis  into  the  right  half  plane,  are  of  radius 
0 < e < e„.  These  semicircular  indentations  are  made  around  all 
points  of  P lying  on  the  imaginary  axis  and  around  all  points  q^ 
of  Q lying  on  the  itu-axis  with  q^  t T.  The  other  "indentations" 
{again  of  radius  e,  0 < e < e0)  are  slits  into  which  encircle  all 

O 

points  of  P and  all  points  q^  of  Q (q^  t T)  which  are  in  £+  so  as 
to  eliminate  these  points  from  the  interior  of  the  contour.  We 
have  also  labeled  these  slits  y-|,  . ..,ur  where  each  maps  an 
appropriate  subtnterval  of  I onto  the  specified  subset  of  (yrK 

O 

The  parellel  lines,  connecting  a pole  in  with  the  semicircular 
portion  of  yr»  are  actually  the  same  line  segment  (slit)  traversed 
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in  opposite  directions.  Note  theft  we  have  indicated  the  usual 
counterclockwise  orientation  to  the  path.  Thus  the  only  points 

O 

encircled  by  yr  are  points  of  T which  are  in  Z+. 


O 
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Let  r denote  the  classical  Nyquist  contour  where 
r:  I -*■  PUt*}  as  indicated  in  Figure  4-b. 


Lemma  2:  Under  the  above  assumptions  on  g and  yr»  h(s)  is  stable 
if  and  only  if  the  path  g»yR  does  not  encircle  "-1". 

A 

Proof:  Since  g(1«*»)  i -1,  -•  <. u <_  •,  there  is  a finite  R such 
that  yR  encircles  all  points  of  T lying  in  jt+.  This  fact  together 
with  g being  analytic  on  X implies  that  the  statements  of  the  Lemma 
are  well  defined.  . 

A A O 

Suppose  h(s)  Is  stable.  Then  h(s)  has  all  poles  in  £_. 

A 

Equivalently  g(s)  j*  -1  for  all  s in  £+.  Thus  yR  does  not  encircle 
any  points  of  T implying  that  yR  is  homotopically  trivial  in  X. 

A 

By  corrollary  to  Lemma  1,  g0yR  is  homotopically  trivial  in  Y. 

A 

Conversely,  suppose  that  g*yR  does  not  encircle  "-1".  Then 

A 

goyR  Is  homotopically  trivial  in  Y.  The  same  corrollary  implies 
yr  does  not  encircle  any  points  of  T.  Thus  all  points  of  X which 

O 

map  to  "-1"  are  in  £_. 


At  this  stage  let  us  compare  the  Information  of  the  Nyquist 

A A 

plot,  g*r,  with  the  "ugly"  Nyquist  plot,  g«Y„  . 

Ro 
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Proof:  Consider  that 


1 ■>  1 -I  r 

a ■ f / _ r /_  i j.  , r 


2*T  / (z-l)‘'dz  * 2*7  / (z-1)  dz  + £ J (z-1)"  dz 

« » i(=i  - 

9oyr  9°r  g®uk 


But  / (z-l)_1dz  - / (g(z)-l)“V(z)dz 

g°uk  uk 

If  uk  encircles  a point  of  Q,  then  (gCz)-!)-1^  (z)  is  analytic 
in  the  region  bounded  by  uk  and  thus  the  integral  approaches  zew 
uniformly  for  arbitrarily  small  e.  Thus  the  Integral  is  zero 
these  points. 

If  yk  encircles  a pole  of  g(s)  then  since  (g(z)-l)“1g’(z)  is 
analytic  in  the  region  bounded  by  uk: 


/ (g(z)-lg '(z)dz  * / d(ln(g(z)-l ) * ln(g(z)-l)  »2irl 

wk  yk  lyk 

for  a suitable  branch  of  the  logarithm.  The  integral  comes  out 

as  negative  2*1  since  pk  was  traversed  in  the  clockwise  direction. 

The  conclusion  of  the  lemma  now  follows. 

At  this  point  let  us  remove  the  restriction  that  g(iw)  f -1 

for  — u <_•.  We  now  give  a proof  of  the  classical  Nyquist 

criterion  using  the  above  concepts. 

A 

Theorem  1:  Let  g(s)  be  as  above  with  the  earlier  restriction 

A 

removed.  Then  h(s)  Is  stable  If  and  only  if  the  Nyquist  plot  of 
g(s)  does  not  pass  through  "-1"  and  encircles  "-1"  exactly  n times 


where  n Is  the  number  of  poles  of  g(s)  in  (L. 

* ' A 

Proof:  Suppose  h(s)  Is  stable,  then  all  poles  of  h(s)  are  in 

A A A 

and  h(»)  < •.  Thus  g(-)  f -1  and  g(s)  f -1  for  all  s in  £+. 
Therefore  via  Lemmas  2 and  3 the  Nyquist  plot  encircles  "-l"  exactly 
n times. 

Conversely  suppose  the  Nyquist  plot  encircles  "-1"  exactly  n 

A 

times  and  does  not  pass  through  "-1".  Thus  g(-)  t -1  which  implies 

A A 

h(»)  < •.  Moreover  since  g»r  encircles  "-1"  n times  and  there  are 

A O A 

n poles  of  g(s)  in  fc+,  we  know  that  g«yR  is  homotopically  trivial. 
Thus  yr  is  homotopically  trivial.  Consequently  there  are  no 

A A 

points  t.  in  £+  such  that  g(t^)  * -1.  Thus  h(s)  is  stable. 


Matrix  Case 


Let  the  entries  of  an  nxn  matrix  G(s)  be  rational  functions 

A 

in  the  complex  variable  s.  Suppose  G(s)  depicts  the  open  loop  gain 
of  the  single  loop  feedback  system  of  Figure  6. 


II 


Figure  5 


x(s)  and  y(s)  are  n vectors  whose  entries  are  also  rational  functions 
of  s which  represent  the  input  and  output  of  the  system  respectively. 


! 1 
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In  this  article,  we  assume  each  entry  of  G(s)  is  bounded  at 

x ■ *.  Thus  G(s)  as  a mapping,  G(*):£  ♦ ?nxn,  is  analytic  on  l except 

at  a finite  number  of  points,  the  poles  of  its  entries. 

For  Figure  6 to  be  well  defined  we  require  that  det  [I+G(s)]  »fc  0 

Thus  there  exists  a closed  loop  convolution  operator,  H,  such  that 

* 

y a H*x.  Moreover  the  Laplace  transform  of  H,  H(s)  satisfies 


H(s) 


G(s)[I+G(s)]“1 


l| 


For  the  system  of  Figure  6 to  be  stable,  H(s)  must  have  all 
its  poles  in  jS_  and  have  all  Its  entries  bounded  at  s * «*. 

A 

Under  the  assumptions  on  G(s),  the  following  factorization  is 
valid: 


G(s)  - N(s)D-1 (s) 

where  N(s)  and  D(s)  are  right  co-prime,  polynomial  matrices  in  s 
with  det[D(s)]£  0.  Moreover  sQ  is  a pole  of  G(s)  If  and  only  if  it 
is  a zero  of  det[D(s)].  (9) 

Oesoer  and  Schulman  (3)  have  shown  that  the  close  loop  operator 
H Is  stable  if  and  only  if  det[N(s)+D(s)]  t-  0 for  s in  and 
det[I+G(")]  £0.  Using  this  fact,  we  state  and  prove  the  following: 
Theorem  2:  H is  stable  if  and  only  if  (1)  the  Nyquist  plot  of 
det[N(s)+D(s)]  does  not  encircle  nor  pass  through  "0",  and  (2) 
det[I+G(-)]  *0. 
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'Proof:  By  hypothesis  we  require  det[I+G(»)]  f 0.  Therefore  we  only 
must  verify  that  det[N(s)+D(s)]  t 0 for  Re(s)  > 0 if  and  only  if  the 
Nyquist  plot  of  det[N(s)+D(s)]  does  not  pass  through  nor  encircle 
”0". 

Now  the  Nyquist  plot  of  det[N(s)+0(s)]  passes  through  "0"  if  and 

only  if  det[N(s)+D(s)]  has  a zero  on  the  Imaginary  axis— i.e.  if  and 

only  if  the  closed  loop  system  has  a pole  on  the  Imaginary  axis. 

Finally  assume  the  Nyquist  plot  of  det[N(s)+D(s)]  does  not  pass 

through  "0".  Observe  that  det[N(s)+D(s)]  is  a polynomial  and  thus 

a rational  function.  As  per  Theorem  1,  appropriately  define  X and  Y 

so  that  X is  a covering  space  of  Y under  the  map  det[N(*)+0(’)]-  The 

above  lemmas  imply  that  the  Nyquist  plot  of  det[N(s)+D(s)]  is  homo^ 

o 

topically  trivial  if  and  only  if  there  exists  a point  Sq  in  £+  such 

that  det[N(s)+D(s)]  * 0.  The  assertion  of  the  theorem  now  follows. 

Observe  that  if  one  assumes  the  open  loop  gain  to  be  stable 

o 

(i.e.  G(s)  has  all  poles  in  £+  ) then  det[I+G(s)]  can  replace 
det[N(s)+D(s)]  in  the  above  theorem.  This  follows  since  for  all  s 
in  fc+,  det[N(s)+D(s)]  = det[I+G(s)]  det[D(s)]  with  det[D(s)]  f 0. 

Thus  in  t+  det[N(s)+D(s)]  has  a zero  if  and  only  if  det[I+G(s)]  has 
a zero. 

Finally,  It  is  worthwhile  to  point  out  the  relationship  between 
the  above  formulated  multivariable  Nyquist  criterion  and  that  formulated 
by  Barman  and  Katznelson!  For  this  purpose  we  let  (iw) ; j»T,  ....  n; 
denote  the  n eigen  values  of  G(iw).  In  general  parameterization  of 
these  function  by  1w  Is  not  uniquely  determined  but  one  can  always 
formulate  such  a function.  Moreover  these  functions  are  piecewise 
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analytic  and  can  be  concatenated  together  in  such  a way  as  to  form  a 

A 

closed  curve  which  Barman  and  Katznelson  term  the  Nyquist  plot  of  G(s). 
Now,  since 


det[I  + G(iw)]  = 11  [1  + A.  (iw)] 

j=l  J 


and  the  degree  of  a product  is  the  sum  of  the  degrees  of  the  individual 
factors  and  also  equals  the  degree  of  the  cancatonation  of  the  factors, 
the  degree  of  the  Barman  and  Katznelson  plot  with  respect  to  "-1" 
coincides  with  the  degree  of  our  plot  with  respect  to  "0".  As  such, 
even  though  the  two  plots  are  different  their  degrees  coincide  and  hence 
either  can  be  used  for  a stability  test. 
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Texas  Tech  University 
Lubbock,  Texas  79409 

Abstract, 

Techniques  for  implementation  of  a Nyqulst  stability  result  for 
a linear  time  invariant  multivariable  feedback  system  are  des- 
cribed. The  approach  is  based  on  continuation  methods  for  com- 
puting the  system's  eigenvalue  loci. 

I.  INTRODUCTION 

The  classical  Nyquist  stability  criterion  for  single-input  single-output, 
linear  time-invariant  feedback  systems  has  only  recently  been  generalized 
to  multivariable  feedback  systems  [1,2].  Stability  theorems  are  expressed 
in  terms  of  the  eigenvalue  loci  of  the  open  loop  transfer  function  G(s)  of 
the  system.  In  particular  if  G(s)  is  stable,  i.e.,  G(s)  has  no  poles  in 
the  right  half  of  the  s-plane  or  on  the  ju>-axis,  then  a linear  time-invari- 
ant multivariable  feedback  system  with  n inputs  and  n outputs  is  stable  if 
and  only  if  its  generalized  Nyquist  plotsfunion  of  eigenvalue  loci)  does 
not  pass  through  or  encircle  the  (-1,  0)  point  [1].  In  order  to  apply  the 
multivariable  Nyquist  criterion,  it  is  thus  necessary  to  compute  the  eigen- 
value loci  as  a function  of  frequency.  For  a given  frequency,  the  eigen- 
values can  be  calculated  by  using  classical  techniques.  Since  the  eigen- 
values are  functions  of  frequency,  normally  one  would  have  to  repeat  the 
entire  computational  procedure  for  each  frequency.  In  the  actual  stability 
analysis,  this  repetition  is  however,  impractical.  Our  approach  to  the 
stability  analysis  of  multivariable  feedback  systems  is  based  on  continu- 
ation methods.  The  basic  idea  of  all  continuation  methods  is  to  convert 
the  solution  of  a parameterized  family  of  algebraic  problems  into  the 
solution  of  a differential  equation.  Then  if  one  can  find  the  solution 
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of  an  initial  problem  by  using  classical  methods  the  solutions  to  the 
other  problems  can  be  obtained  by  integrating  the  associated  differen- 
tial equation  with  the  initial  solution  as  an  initial  condition. 

II.  EIGENVECTOR  APPROACH 

Our  first  method  is  based  on  the  approach  described  by  Faddeev  and 
Fadeeva  [3]  and  Van  Ness  et.  al.  [4].  A differential  equation  is  written 
with  the  eigenvalues  as  dependent  variables  and  the  frequency  as  variable 
parameter.  We  then  compute  a set  of  initial  eigenvalues  by  classical 
analysis  techniques  and  integrate  the  resulted  differential  equation  to 
obtain  the  required  eigenvalues  for  each  frequency.  The  eigenvalues  x^(o») 
of  G(joi)  and  their  complex  conjugates  x.(u)  satisfy 

G(jw)X.U)  = X^oOX-U)  i=l  ,2, . . . ,n  (1) 

and 

G (ju>)V.(u)  = X.UJV.U)  i-l,2,...,n  (2) 

where  X.(w)  and  V.(u)  are  the  corresponding  eigenvectors  of  x.(w)  and 
x.j U)  respectively,  and  G (j<j)  is  the  complex  conjugate  transpose  matrix 
of  G( jtu) . 

We  differentiate  (1)  with  respect  to  u to  yield 


xi  V 


“j  i “ 1 » 2>  • • • » n • 


da>  <X . , V •>  ’ ' B#  - 

i i 

The  differential  equations  involving  and  are  obtained  as 


r-jii  1=1,2 

dVi  - l 8.-.V.,  i = 1,  2,  ....  n. 
d ~ ~ j-1  10  J 


where 


*dS  V V i i j. 

“11-°-  “1J-  (»rVJ7dpj» 


3n  * °-  6ij 


<V  Y 


i t j. 


Starting  with  a set  of  predetermined  initial  conditions  x.j(0)  = xiQ, 
X;-(0)  = X.Q  and  V.(0)  = V^Q  for  i =1,  2,  ....  n,  we  integrate  (3),  (4) 
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ar|d  (5)  to  obtain  the  required  eignevalues  for  each  frequency.  The  eigen- 
value loci  are  computed  in  a continuous  manner  by  numerical  integration. 

III.  JACOBIAN  METHOO 

For  an  nth  order  system,  the  above  algorithm  requires  the  numerical  inte- 
gration of  a set  of  3n  equations  and  the  computation  of  two  sets  of  un- 
wanted variables— namely  the  eigenvectors  X.  and  V^.  These  disadvantage, 
can  easily  be  avoided  if  the  characteristic  equation  for  the  multivariable 
feedback  system  can  be  predetermined.  A much  simpler  method  can  be  formu- 
lated based  on  the  approach  for  finding  multiple  solutions  for  a nonlinear 
equation  developed  by  Chao  et.  al.  [5]. 

Let  the  characteristic  equation  of  G(ju)  be  given  by  an  nth  order  polynomial 
in  eigenvalue  x with  complex  coefficients 

f|>U)]  - I Al  - G(ju>)  | * 0.  (8)  1 

Instead  of  solving  (8)  directly  for  each  frequency,  we  consider  two  simul- 
taneous differential  equations  of  the  form 


jj£-  -f(t)  f(0)  = f[A(u0)]  = 0 

3t  " * 1 “(0)  3 V 

Assuming  the  nonsingularity  of  the  Jacobian  Matrix 
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in  the  x-«u  space  the  algorithm  (9)  reduces  to 
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It  is  seen  from  the  solution  of  (9) 
f(t)  - 0e“t  = 0 
u * + t. 


(9) 
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that  for  any  admissible  pair  of  uQ  and  ,\(o>Q)  satisfying  (8),  the  corres- 
ponding trajectory  will  remain  on  the  solution  curve  f=0  as  u changes.  The 
+ or  - sign  is  chosen  depending  on  whether  one  would  like  to  increase  or 
decrease  w.  Equation  (11)  may  now  be  solved  by  any  numerical  integration 
techniques  and  the  eigenvalue  loci  can  be  traced  automatically  by  inte- 
grating only  a second  order  differential  system. 

IV.  EXAMPLE 

To  illustrate  the  approaches  presented,  consider  a linear  time-invariant, 
multivariable  feedback  system  with  open  loop  transfer  function  character- 
ized by 


G(s)» 


4 

s+2 

s+T 


k 

s+2 

4 


03) 


for  which  the  characteristic  equation  is  given  by 
f[x(«)]  = x2-8x  + (--a  + j25w_  )aQ.  (14) 

1 + a.  1 + </ 


The  generalized  Nyquist  plots  shown  in  the  accompanied  figure  for  the 
cases  where  k*16  and  36  are  obtained  by  applying  the  eigenvector  approach 
where  as  in  the  critical  case,  k=25,  the  Jacobian  method  has  been  used. 
In  all  three  cases,  the  equations  are  integrated  using  Euler's  method 
with  a step  size  of  0.01.  It  is  seen  from  the  figure  that  the  system  is 
stable  for  k<25  since  the  generalized  Nyquist  plots  do  not  encircle  -1 
point. 
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